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Abstract. The two-phase free boundary value problem for the isothermal 

Navier-Stokes system is studied for general bounded geometries in absence of 

phase transitions, external forces and boundary contacts. It is shown that the 

problem is well-posed in an Lp-setting, and that it generates a local semifiow on 

P^ i the induced phase manifold. If the phases are connected, the set of equilibria of 

^M ' the system forms a (ra+l)-dimensional manifold, each equilibrium is stable, and 

^J it is shown that global solutions which do not develop singularities converge 

r^ , to an equilibrium as time goes to infinity. The latter is proved by means of 

the energy functional combined with the generalized principle of linearized 
stability. 
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O . 1- Introduction 

In this paper we consider a free boundary problem that describes the motion 
of two isothermal, viscous, incompressible Newtonian fluids in R'^. The fluids are 
^^ . separated by an interface that is unknown and has to be determined as part of the 

?H ' problem. 

^ ■ More precisely, we consider two fluids that fill a region il C M^ . Let Fg C fJ be a 

given surface which bounds the region fii(O) occupied by a viscous incompressible 
fluid, fluidi, the dispersed phase and let ^2(0) be the complement of the closure of 
Oi(0) in Q, corresponding to the region occupied by a second incompressible viscous 
fluid, fluid2, the continuous phase. Note that the dispersed phase is assumed not 
to be in contact with the boundary dfl of fl. We assume that the two fluids 
are immiscible, and that no phase transitions occur. The velocity of the fluids is 
denoted by u{t,x), and the pressure field by 7r(i,x). 

Let r(t) denote the position of Fq at time t. Thus, F(i) is a sharp interface which 
separates the fluids occupying the regions ^\{t) and 1^2 (i), respectively. We denote 
the normal field on F(i), pointing from J7i(i) into ^2{t), by j/r(i, •)• Moreover, 
Vr(i, •) and Hr{t, ■) mean the normal velocity and the curvature of F(t) with respect 
to i^rit, ■), respectively. Here the curvature Hr := — divrJ^r is negative when fli(t) 
is convex in a neighborhood of a; G F(t), in particular the curvature of a sphere 
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Sr{xo) is —{n — l)/R. The motion of the fluids is governed by the fohowing system 
of equations for i = 1,2. 

Pi(9(U + (u|V)u) — /iiAu + Vtt = in ^i{t), 

divu = in ^i{t), 

— |S'(m, 7r)t'r] = crHrt^r on T{t), 

M=0 on r(t), 

u = on dVt, 

Vr = iu\iyr) on T{t), 
u{0) — uq in fi\ro, 

r(o) = To . 

Here, 5* is the stress tensor defined by 

S = S{u,tt) ^ fii{Vu+[Vu]^) -nl ^2iJ,,E -7:1 in n,{t), 

and 

Ul{t,x)= hm ((j){t,x + hiyr{x)) -(l)(t,x -hi^rix))), xeT{t) 

denotes the jump of the quantity (/>, defined on the respective domains fli{t), across 
the interface r{t). 

Given are the initial position Fq of the interface and the initial velocity uq : r2\ro — >■ 
M.^. The unknowns are the velocity field u{t, ■) : 51\r(i) — > M'^, the pressure field 
TT{t, •) : n\T{t) -)■ K, and the free boundary T{t). 

The constants pi > and pi > denote the densities and the viscosities of the 
respective fluids, and the constant a > stands for the surface tension. In the 
sequel we drop the index i since there is no danger of confusion; however, we keep 
in mind that /x and p have jumps across the interface, in general. 

System (II. 1[) comprises the two-phase Navier-Stokes equations with surface ten- 
sion. The corresponding one-phase problem is obtained by setting P2 = P-2 = and 
discarding J72. Here we concentrate the discussion on the two-phase problem. 

There are several papers in the literature dealing with problem (jl.ip ; cf. [5j |6j 
[71 [HI m [13 Uni [23 ■ All of them employ Lagrangian coordinates to obtain local 
well-posedness. This way it seems difficult to establish smoothing of the unknown 
interface, and this method is hardly useful in case phase transitions have to be taken 
into account. Here we employ a different approach, namely the Direct Mapping 
Method via a Hanzawa transform, which has been quite efficient in the study of 
Stefan problems, i.e. phase transitions involving temperature, only. 

In a recent paper [19] we have shown that problem (jl.ip is locally well-posed in 
an Lp-setting provided H. = K" and the initial interface Fq is sufficiently close to 
a flat configuration. In addition, the interface as well as the solution are proved 
to become instantaneously real analytic. This result is based on a careful analysis 
of the underlying linear problem. Building on the latter results we show in this 
paper local well-posedness for arbitrary bounded geometries as described above. 
This induces a local semifiow on a well-defined nonlinear phase manifold. 

It is known that the set £ of equilibria of the system are zero velocities, constant 
pressures in the components of the phases and the dispersed phase is a union 
of disjoint balls. Concentrating on the case of connected phases, we prove that 
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equilibria are stable and any solution starting in a neighbourhood of such a steady 
state exists globally and converges to another equilibrium. 

The energy of the system serves as a strict Ljapunov functional, hence the limit 
sets of the solutions are contained in the set of equilibria £. Combining these 
results we show that any solution which does not develop singularities converges to 
an equilibrium in the topology of the phase manifold. 

2. Transformation to a Fixed Domain 

Let 51 C M" be a bounded domain with boundary dVt of class C^ , and suppose 
r C ri is a hypersurface of class C^, i.e. a C^-manifold which is the boundary of 
a bounded domain J7i C fl; we then set fl2 = fl\fli. Note that fl2 is connected, 
but ill maybe disconnected, however, it consists of finitely many components only, 
since dili = F by assumption is a manifold, at least of class C^. Recall that the 
second order bundle of F is given by 

Here Vr denotes the surface gradient on F. Recall also the Haussdorff distance dn 
between the two closed subsets A,Bc M™, defined by 

duiA^ B) := max{sup dist(a, B), sup dist(6, A)}. 

aeA beB 

Then we may approximate F by a real analytic hypersurface S, in the sense that 
the Haussdorff distance of the second order bundles of F and E is as small as we 
want. More precisely, for each 77 > there is a real analytic closed hypersurface S 
such that d/f (TV^S, A/'^F) < r/. If 77 > is small enough, then E bounds a domain 
Gi with Gi C ri, and we set G2 = fl\Gi. 

It is well known that such a hypersurface S admits a tubular neighbourhood, 
which means that there is a > such that the map 

A : E X (-a, a) ^K", A(p,r) := p + riy^{p), p G E, |r| < a 

is a diffcomorphism from E x (—a, a) onto TZ(A). The inverse 

A-^ :7^(A) K^ E X {-a,a) 

of this map is conveniently decomposed as 

A-^{x) = (n(x), d{x)), X e 7^(A). 

Here H(a:;) means the orthogonal projection of x to E and d{x) the signed distance 
from X to E; so \d{x)\ = dist(a::, E) and d{x) < iff a; G Gi. In particular we have 
7^(A) = {x e K" : dist(x, E) < a}. 

Note that on the one hand, a is determined by the curvatures of E, i.e. we must 
have 

< a < min{l/|Kj(p)| : j = 1, . . . , n — 1, p G E}, 

where Kj{p) mean the principal curvatures of E at p G E. But on the other hand, a 
is also connected to the topology of E, which can be expressed as follows. Since E 
is a compact manifold of dimension n — 1 it satisfies the ball condition, which means 
that there is a radius r^ > such that for each point p G E there are Xj G Gj, 
j — 1,2, such that Br-^{xj) C Gj, and Br-^{xj) n E = {p}. Choosing r^ maximal, 
we then must also have a < r-^.- 
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Setting r(0) = Fq, wc may use the map A to parametrize the unknown free 
boundary T{t) over S by means of a height function h via 

r{t) = nip^p + h{t,p)i^^ip), pes), t>o, 

for small t > 0, at least. Extend this diffeomorphism to all of f2 by means of 

Ohit, x)=x + x{d{x)/a)h{t, n(x))i.s(n(x)) =: x + 0,,(i, x). 

Here x denotes a suitable cut-off function; more precisely, x G 2^(1^) j < x < 1, 
x(r') = 1 for \r\ < 1/3, and x(^) — for l^j > 2/3. This way il\r(i) is transformed 
to the fixed domain il\S. Note that &h(t,x) = x for \d{x)\ > 2a/3, and 

Q-\t,x)^x-h{t,U{x))i^^{Il{x)) for |d(a;)| < a/3, 

in particular, 

Q^ {t,x) = X — h{t,x)i's{x) for a; G S. 
Now we define the transformed quantities 

u{t,x) = u{t,Qhit,x)), 
Tt{t,x)^Tr{t,ehit,x)), i > 0, a; e f7\E, 

the pull backs of u and tt. This gives the following problem for u, tt, h. 

pdtu - fiA{h)u + g{h)Tf = n{u, h) in f]\E, 

{Q{h)\u) = inf7\S, 

^l^l{[g{h)u] + [g{h)uY) - v:}vT{h) = aHrihyrih) on S, 

|u] =0 on S, (2.1) 

u — on 9r2, 

9tft- — (u|t^s) = ^{u\a{h)) on S, 

u(0) = uq, in ri\S, /i(0) = /iQ on S. 

Here ^(/i), ^(/i) and Hr{h) denote the transformed Laplacian, gradient and cur- 
vature, respectively. More precisely, we have 






and 



and 



e',^i + el e'^-' ^ I - [I + e',] 
[Vtt] o e,, = gih)n 

=: (I -Mi{h))VTt 

[divu]oe^ = {g{h)\u) 

= mn'f ° ejv|s) = {v\u) {e';![i + e',]-^v\u) 

= {{I - Mi{h))^\u) 
[Au] 06/1 = ^(ft,)u 

= [(Qh ')'(0h ')'"" ° Q'^] : v^s + ([Ae-1] o e;,|v)K 

= Au - M4(/i) : V^u - M2{h)Vu 
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with 

-M2{h)\/u := ([Ae,;i] o e,,|v)fi 

M,ih) : S7^u := [2symi9';^[I + 9',]-^) - [I + e',^r'e',^e'^[I + 6',^^] : V'u. 
Note that 

[dtu] o e,, = dtu - u'[{dte^^) o e,,] = a^^ - s'e;riateft 
= a^M - u'[i + e'l^r^e'i^dtOh =■. dtU - M3ih)Vu, 

hence 

R{u, h) = -p{u ■ g{h)u) + M3{h)\/u. 

With the curvature tensor L^: and the surface gradient Vs we have 

Mh) - /3(M(f^s - a(M), a(/i) = Moih)V^h, 

Mo(h) = {I- hLj:)-\ /3{h) = (1 + |a(/i)|2)-i/2^ 
and 

Employing this notation, we have 

0'hit, x) = i^E(n(x)) Mo{d{x))V^hit, U{x)) 

- h{t,nix))L^(fl{x))Mo{d{x))V^ for |d(a;)| < a/3, 
ehit,x) =0 for |d(a;)| > 2a/3, 



and 



9h{t,x) - ix'(c?W/a)/i(i,n(x))i.E(n(x)) ® i.s(n(a;)) 
a 



+ xidix)/a)i^^iU{x)) (E) Mo{d{x))Vj:h{t,U{x)) 
- x{d{x)/a)h{t, n(x))Ls (n(x))Mo {d{x))V^ 
for a/3< |d(x)| < 2a/3, 

where Ps = I — v^®vy^ denotes the projection onto the tangent space of S. Thus, 
[/ + 9'f^ is boundedly invertible, if |ft,|oo and |Vs^|oo are sufficiently small. The 
curvature Hr{h) becomes 

Hrih) ^ Pih){tT[Moih){L^ + VMh))] - l3^{h){Mo{h)aih)\[Vj:a{h)]a{h))}, 

a differential expression involving second order derivatives of h only linearly. Its 
linearization is given by 

i?f(0) = trL| + As. 
Here As denotes the Laplace-Beltrami operator on S. 

It is convenient to decompose the stress boundary condition into tangential and 
normal parts. Multiplying the stress interface condition with ly-s/fi we obtain 

H - aHrih) - {M[g(h)u] + [g{h)u]'')j{iyj: - Mo{h)WshWs), 

for the normal part of the stress boundary condition, and 

- P^M[G{h)u] + [g{h)u]')j{iy^ - Mo{h)W^h) 

= {M[g{h)u] + [g(/i)w]T)](j.s - Mo(/i)Vs/i)|^s) Mo(h)V^h, 

for the tangential part. Note that the latter neither contains the pressure jump nor 
the curvature! 
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We rewrite this problem in quasilinear form, dropping the bars and collecting its 
principal linear part on the left hand side. 

pdtu - fiAu + Vtt = F{h, u)Vu + AUih) : V^u + Mi(/i)V7r in n\T,, 

divu^ Mi{h) -.Vu in f^\S, 

P^l-fi(Vu + Wu'^)jiyj: = Gr(h)Vu, (2.2) 

-{l-^i{Vu + Vu^)jiy^\iy^) + H - aiJf (0)/i = G^(/i)Vm + G-,{h) on E, 

|u] = on S, 

u — on dfl, 

dth - {u\v^) = (Mo(/i)Vs/i|u) on E, 

u(0) = uo in i^\S, ft,(0) = /iQ on S. 

The right-hand sides in this problem are either lower order terms or are of the same 
order appearing on the left, but carrying a factor h or Ve/i, which are small by 
construction. In fact, since Fq is approximated by S in the second order bundle 
we have smallness of ft-o, Ve/iq, and even of V|;/io, uniformly on E. All terms on 
the right-hand side are at least quadratic. More precisely, besides the Mj{h) which 
have been introduced before, the nonlinearities have the following form: 

F{h,u)SJu = - (u|Vu) + [Mi{K) + A'hih) + M:i{h)]Vu, 
G,{h)yu = - iMlVu] + [Vu]'^)]Moih)Vsh\iys) 

- {ln([Ah{h)Wu] + [Mi(/i)V?/]T)](j.E - Mo{h)W^h)\iy^), 

Gr{h)\/u =(|m([(/ - Afi(/i))Vu] + [(/ - Mi(/i))Vm]T)](i.s - Mo(/i)VsMl^s)- 
■Moih)V^h 

- P^Miil - Ah{h))Vu] + [{I - Ahih))Vu]^)lAIoih)Vj:h 

- P^lti{[Ah{h)yu] + [Mi(/i)Vu]T)]i/E, 
G-,{h) ^ (j{Hr{h) - H{.{0)h). 

The idea of our approach can be described as follows. We consider the transformed 
problem (12.21) . Based on maximal Lp-regularity of the linear problem given by the 
left hand side of (j2.2p , we employ the contraction mapping principle to obtain local 
well-posedness of the nonlinear problem. The solutions of the transformed problem 
will belong to the following class: 

u e i/pi(J;ip(r!)") nLp(J;i72(r)\E)"), t: G Lp(J; ffpi(f]\E)), 

h G H^p2-i/2p( J. ^^(^)) ^ fji(^j. w^-^^'^ij:)) n Lp( J; Wp3-i/''(S)). 

This program will be carried out in the next sections. 
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3. The Linearized Problem 
We consider now the inhomogeneous linear problem 

pdtu — fiAu + Vtt ~ pf in r2\I], 

div u = fd in 51\I], 

l-piVu + [Vw]"^) + iTTpE, - ct(Ae/i)i's = g on E, 

M = lis on E, (3.1) 

u ^ ui, on dfl, 

dth - {u\iy^) + (6|Vs/i) = gh on E, 

ii(0) = uq in r2\E, /i(0) = /iQ on E 

on a finite time-interval J = [0, a]. We choose the same regularity classes for u and 
TT as before, i.e. 



and 

Then 

and 



ueZ^:= H^{J; Lp{nr) n Lp{J; ij2(f7\S)"), 



ub G VFp^"^/'p(J; ip(9n)") n Lp{j; w^-'^p{dnr). 

Therefore the equation for the height function h lives in the trace space for the 
components of u, i.e. 

5, G y„° := W^''/'^{J; ip(E)) n ip( J; M^^i/^(E)), 

hence the natural space for h is given by 

Here the last space comes from the curvature term in the stress boundary condition, 
which induces an additional order in spatial regularity. Assuming that g belongs 
to the trace space of Vu, i.e. 

g G rj := w;,/^-'/^p{j;Lp{j:r) n Lj,{j;w;,-'^p{j:r), 

we have the additional regularity [tt] G Y^ for the pressure jump across the interface 
E. The function b G Y^ is given; we will choose b appropriately in Section |4l 

There is another hidden regularity which comes from the divergence equation. 
To identify it, let (j> G HL{n). An integration by parts yields 



{u\V4>)n = -{d[vu\(j))n + (u ■ J^aa|0)9n - d" • '^T,j\4>)^ 
= -{fd\(t))n + {ub ■ ydn\(i))dn - ("e • J^e|'/')s- 

Set H~'^{Vl) = {Hp,{fl))* and define the functional {fd,Ub ■ Vdn,us ■ i^t.) G Hp^{fl) 
by means of 

{{fd,Ub ■ VQQ,,UT. ■ VT.)\(i)) := -{fd\(t>)n + [ub ■ vdn\(i>)dn - (us • VT.\(i>)T.- 

Then we have 

{{fd,ub ■ vdn,UT, ■ VT,)\4>) = (u|V(/))a. 
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Since u G H^{J; Lp{n)"-) this implies (/d,U6 • Vdn^u^ ■ v^.) G H^{J;Hp^{n)). Ob- 
serve that this condition contains the compatibility condition 



fddx ^ Ub ■ i/on ddVL - I ut, ■ i^t, dY., 
n Jon Jt. 

which appears choosing 0=1. 

In the particular case fd — Owe have {fd,Ub ■ vdn,UY, ■ v^) G Hp{J; Hp^{i})) if 

and only if Ub ■ van e Hl{J- ^-^'"^[8^)) and u^ ■ y^ ^ H^{J-, T4V'^''(S)). 

The main theorem of this section states that problem p.ip admits maximal 
regularity, in particular, it defines an isomorphism between the solution space and 
the space of data. 

Theorem 3.1. Let p > n + 2, fl C M" a bounded domain with dfl G C^ , S C fi 
a closed hypersurface of class C'^ and pj , fij, a he positive constants, j — 1,2; set 
J = [0,a], and suppose 

b G W^-^/^P{J; Lp{j:)r n Lp{J; M/p2-i/P(s))». 

Then the two-phase Stokes problem iS.l]) admits a unique solution (w, tt, h) with 
regularity 

u G H\j-Lp{nr)r\Lp{J-Hl{n\i:Y), vr G Lp{J;H^{n\j:)), 

h G W^-^/'P{J; Lp(S)) n H;,{J; W-,2-i/p(s)) n Lp(J; W^-^/p{^)), 

if and only if the data (uq, hQ,Ub,u-s, f, fd, g, gh) satisfy the following regularity and 
compatibility conditions: 

(a) f G Lp{J X nr, uo G W^'^^^inXJ:)", and uqi,,, = Ub,^^„; 

(b) fd e Lp{J; Hl,{n\Y)), and dW uo = fd^^,; 

(c) Ub G I<^'/'P( J; Lpi^d^Y) n ip( J; Wp-^'^^id^r), 

(d) «E G Wl-'/'^iJ; Lpi^r) n Lp{J; W^p'-^/^(I])«); 

(e) {fd,Ub ■ VdQ,,u^ ■ vt) e H^{J;H~^{n)); 

(f) g G wl'^-"^^{J-Lp{i:)rnLp{J-W^-"^{'i:)r; 

(g) IM = «s|.^„, and Ps[/i(VMo + [Vu^Y)} = P^g^^,; 

(h) ho G W^-^^^i^), and gn G W^^'^'^J-Lpi^)) H Lp(J; ^^'-^/'^(S)). 

The solution map {uQ^ho,Ub^u^^ f, fd, g, gh,b) f-> (w, tt, [tt], /i) is continuous be- 
tween the corresponding spaces. 

The proof will be carried out in the following subsections. 
In general the pressure w has no more regularity as stated in Theorem l3.1l However, 
there are situations where tt enjoys extra time-regularity, as stated in the following 



Corollary 3.2. Assume in addition to the hypotheses of Theorem \3.1\ that 
uq = ho = fd = 0, div f = in r2\I], 

Ub ■ von — on dil, u^ ■ v^ — on E, 
and 

[(/I^e)1=0o«S, ifluon)^ Gondii. 

Then n G oH^{J\Lp{VL)), for each a G (0, 1/2 - l/2p). 
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Proof. Let g E Lpi{Vt) be given and solve the problem 

A0 — pg in 51\I], 
101 = on E, 
Ip-^d^cj,} =0 on E, (3.2) 

d^(j) — on dfl, 
by Theorem 18. II Since (/|V0) = (u|V0) = we obtain by integration by parts 

{n\g)n - (-\aA =- f [-^..^IdE - f — |V<^) 
\P Jn Jt. P \ P Ja 

l^'^s'^jv. (P 



^s P \P Jn 

Jn P JT. P P 

Since Vu G oi^y'(J; ip(fi)"^") and H,9feUz G oW^p'/'-'/'^(J; ^p(S)), applying 
i9f to this identity, we obtain the estimate 

for each a G (0, 1/2 - l/2p), hence n G o^p (J; ip(^))- □ 

It is convenient to reduce the problem to the case 

uo = ho = f = fd = u^ ■ v^ = Ub ■ van = 0. 

This can be achieved as follows. Suppose (w, tt, h) is a solution of p.ip . We intro- 
duce a further dummy variable q := |7r|; note that q £ Zq := Y^ . We decompose 
u = u^, + ui, TT = TT^, + TTi, q = q^ + qi, h = h:t + hi where 

h^t) = [2e-(^-^-)'''* -e-2(^-^-)'''>o+ 

[g-(/-A,)t _ g-2(7-A,)t](^ _ As)-H(«o|^s) - (&|Vs/i) + .9^(0)}, t > 0. 

The function ft,, belongs to Z^ and satisfies /i,(0) = ft-o and 9t/i,(0) = (uo|i^s) — 
(fe|Vs/i) + .g/i(0). Then hi has initial value zero, and also dthi{0) = 0. We set 
q,(i) = e'^^*qo where 

go := ([m(Vuo + [yuo]'^)P^\iy^) + aA^ho + {g{0)\i^^) 
is determined by the data, and wc define vr, as the solution of 

Att, = in n\j:, 

9iy7r, =0 on dV,, 

I(9i.j,^*l = 0, [tt,] = g, on E. 

Note that g, G yj and tt, G ^tt, by Theorem 18.51 The function u, G 2^^ is defined 
as the solution of the parabolic problem 

pdtu — p,Au = — Vtt, + pf in ri\E, 

u — Uf, on (917, 

|-M(Vu+[Vu]"^)]i^s=5-9*'^s+CT(Asft,)i^s on E, (3.3) 

|u| == Us on E, 

u(0) = Mo, in r2\E, 
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which is uniquely solvable since the appropriate Lopatinskii-Shapiro conditions are 
satisfied; see [12]. Thus we may assume w.l.o.g. uo — ho = f = and that the 
time traces of fd, g and g^ are zero at time zero. Finally, to remove fd, we solve 
the transmission problem 

A^^fd in f7\I], 

|p^]=0 on S, 

[9,,V']=0 on S, 

^!'anV' = on dVl, 

according to Theorem 18.61 where fd := fd — div u*. Since d^ G C^ , the solution 
satisfies Vi/' € Z^. Then setting U2 = ui — Vij) and 1^2 = tti+ pdtijj — fJ-A.^^ /12 = hi, 
we see that we may assume also fd = uy.-vt. = Ub- ven = 0, the only no n- vanishing 
data which remain are g,gh,UY.,Uf,\ note that the time traces at t = of these 
functions are zero, and Uf, ■ vqq, = on 957 and us • i^s = on S. 

3.1. Flat Interface. In this subsection we consider the linearized problem for a 
flat interface. 

pdtu- ^iAu + Vtt = pf in M", 

div u^ fd in 

-{pdyvj - ln'^xwj = gy on 

~2lpdyWJ + lnj-aAh = gy, on R"-\ (3.4) 

M = us on M"-\ 

dth-w + {b\Vh):^gh on R"-\ 

u(0) = Uo, h{0) = ho in R", on R"-\ 

Here we have identified W"--^ == M"-i x {0} and M" = M"\R"-i. It is convenient to 
split u = {v,w), f = {fv,fw), g — {gv^g-w) into tangential and normal components. 
The following result, which is implied by j20[ Theorem 3.1], states that prob- 
lem p.4p admits maximal regularity, in particular defines an isomorphism be- 
tween the solution space Z := Zu x Z^^ x Zq x Zh and the product-space of data 
(uq, ho, Us, /, fd, g, gh, b) which we denote for short by Y. 

Proposition 3.3. Let p > n + 2 be fixed, and assume that pj, pj, a are positive 
constants for j = 1, 2, and let J = [0, a]. Suppose 



pn 



Dfi— 1 



and set b — bo + bi. Then the Stokes problem with flat boundary p.4p admits a 

unique solution 

(u, TT, h) with regularity 

u e H^{J;Lp{R^r)nLp{J;H^{WT), ^ e Lp{J;H^{W')), 

h e W^-^/^P{J; Lp(R''-^)) n H^{J;W^-^/P{R''-^)) n LpiJ;W^"^/PiW'~^)) 

if and only if the data {f,fd,g,gh,uo,ho,u^) satisfy the following regularity and 
compatibility conditions: 
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(a) / e Lp(JxR")", us e VKp~'/'^(J;Lp(K"-i)")nip(J;W^p"'/^(M«-i)"), 

(b) fd e Lp{J;H^{R^)), {fd,u^ ■ z^s) G Hl{J;Hp\W^)), 

(d) 9h e T<-^/'P(J;Lp(R"-i)) nLp(J; Wf-^/''(M"-i)), 

(e) Mo e VKp"^/P(]R")", /io e Wp"^/^(M"-i), 



(f) divuo = /d|t=o «"• 



id |uo] = us|t=o o?i 



(g) -{^idyVol - IplVxWqI =gv\t=o on K" ^ 

T/ie solution map [{uo,ho,us, f, fd, 9, QhTb) H> (u, tt, |7r], ft,)] is continuous between 
the corresponding spaces. 

3.2. Bent Interfaces. Next we consider the case of a bent interface. By this we 
mean a situation where the interface S is given as a graph of a function ((> : M"^^ — > 



of class BC ; thus E — {{x, 0(x)) : x G 



}. The normal vs is then given by 



iy^{x) = I3{x) 



'Voo(j){x) 

1 



/3(x) = l/Vl + |V,(/.(a;) 



and the Laplace-Beltrami operator for such a surface with 

h{t,x) = h{t,{x,(j){x))) 

reads as 

As/i ^Ah^ /3^{W^hW(l)\V4>) - I3^[A(I) - l3^(V^(l)W(l)\V(t>)]{W(l)\Vh). 

We may assume by the reduction explained above that uq = Hq = f = fd = 
uj: ■ vj: = 0. Set 

u{t, x, y) ^ u{t, x,y + <j>{x)), 7f (t, x, y) = 7r(t, x,y + (l){x)), 

for t E J = [0,a], X ^ M"^^ , y 7^ 0, and observe 

Vu = Vu — V(j) ® dyii. 

Then we obtain for the new variables (u, tt, h) the following problem. For conve- 
nience we drop the bars, and split u ~ {v, w) and g = {gy, g^) as before. 

pi9tu - /xAu -h Vtt = yuBi (u, tt) in K", 

divu = B2U in 

-{fidyvj - IfiVxivj = gy + Bz{u, [■K\,h) on 

-2l^i^yw} + [tt] - crA^/i = [g^//3] + S4(w, h) on 

|u] = ue on 

dth — w + (6|Vft) = gh + B^u + Be/i on 

u(0) = 0, /i(0) = in 



r, 



(3.5) 
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Here we have set 

Bi{u,7r) = \V(t)\^d^u-2{V(j)\V^dyu) + {V(j))dyTT - {A(l})dyU 

B2U = {y<p\dyu), 

B^iu, Inlh) = -l^iiV^v + [V,v]')jV^ - l^idyVJ\V4^f 

+ {-Mdyv\Wcl>)j + H - Ifidywj - aAs/i}V(/) 
Bi{u,h) = -iMdyV + V,zi;)l|V0) - l^idyWJlV^l'' + a{A^h - Ah) 

B,u = (/3 - l)w - p{V(t>\v) - - ^ll^p ^ - Pi^M 
Beh = /32[(fe|V<^) - (6|e„)|V<^p](V0|V/i). 



Now suppose (u, tt, h) belongs to the maximal regularity class. We estimate the 
perturbations Bj as follows: 



+ \A<J)\lJVu\l^, 
|S2u|l,(H1) < |V<^|l^ NMl, + (|V20|l^ + \V^\lJ\Vu\l„ 
l^*^2u|^^(^-i) < |V<?!)9tu|ip(ip) < \S7<l)\Ljdtu\L^, 

\B^{u,M,h)\w^^L,)<C\V<P\L^{l + \V<P\Lj[\Vu\w.^L,) + \V'h\w.^Lj 

\B4{u,h)\wf(^L,)<C\\/c^\L^{l + \\V(f>\\Lj[\Vu\w.^L,) + N^h\wf(^Lj 



Here C denotes a constant only depending on the parameters /i and ct, and we have 
set s = 1/2 — l/2p. For the estimations in Lp{J; Wp iW^^^)) we observe that 



]^^n~l+sp■ 



H\wi^\i^K + ms,„ m.p= I I mx + h)-Mx)\p-^'' 

'\h\<l ' 



defines a norm on Wp{M." ^). This implies 



\atp\w^ < \a\L^\ip\wi +Cs,p|Va|i^|V'|Lp, 
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for a G Wl^iW^^^), with some constant Cs,p which only depends on s e (0,1), 
p G (1, cxd) and on n. With this observation we have 

+ |V20|L^[|Vw|i, + IWIl, + \Vh\L,(Hl)]} 
\B^{u,h)\LAW-^)<C{l + \V<j}\Lj{\V<j,\LjVu\L^(W-s) 

+ C{\V^^\l^ + |V20|D|V0U^|V/i|i,}. 

Here C denotes a constant only depending on /i, a, p, and 2s — 1 — 1/p. To estimate 
Bgh we note that Y^ is a Banach algebra since p > n + 2. This yields 

|B6/i|y^o < C\V(b\L^{\bo\ + \bi\Yo)\yh\YO < C|V0|l^(|6o| + |6l |y„o ) |/l|z, . 

To solve the problem p.Sp . let z = (u,7r, |7r], /i) G O'Z^, where qZ means the 
solution space with zero time trace at i = 0, F := (0, 0, 5^,^10 //3, "s, 5/0 £ o^^, the 
space of data with zero time trace, and let i? : o-^ — >■ o^ defined by 

Bz = {Bi{u, tt), B2U, B^iu, H, /i), Bi{u, h), 0, B^u + B^h). 

Denoting the isomorphism from 0^ to o^ defined by the left hand side of p.Sp by 
L we may rewrite problem p.Sp in abstract form as 

Lz^Bz + F. (3.6) 

The above estimates for the components of B imply 

\Bz\y < C\V(ly\Ljz\z + M[\u\l^(^hI) + \M\l, + \Vh\w^(^L,)nL,(Hl)l 

with some constants C > depending only on the parameters and Af > 0, which 
depends also on | V0|st/c^ ■ Let 77 > be given and suppose | V^JLo^ < "H- By means 
of an interpolation argument we find a constant 7 > 0, depending only on p such 
that there is a constant Af (77) > such that 

\Bz\Y<C[2r] + a?M{ri)]\z\z, z^^Z. 

Choosing first 77 > and then a > small enough, we can solve (13.61) by a Neumann 
series argument for J — [0, a]. 

Since problem (J3.5I) is time-invariant, we may repeat these arguments finitely 
many times, including the reduction procedure, to solve p.Sp for J = [0,a], where 
now a > is arbitrary. 



3.3. General Bounded Geometries. Here we use the method of localization. 
By assumption, d^ is of class C^ and E will even be real analytic, so in particular 
of class C^. Therefore we may cover E by A^ balls Br/2{xj) with radius r > 
and centers Xj G E such that E n Br{xj) can be parameterized over the tangent 
space Tx-T, by a function 9j G C^ such that \Vdj\L^ < r], with 77 > defined as 
in the previous subsection. We extend these functions 9j to all of Tj; . E retaining 
the bound on V^j. This way we have created N bent half-spaces Ej to which the 



(3.7) 
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result proved in the previous subsection applies. We also suppose that Br{xj) C il 
for each j. Set U := il\ Ui=i Br/2{xj) and Uj — Br{xj), j = 1, . . . , N. The open 
set U consists of one component Uq characterized by dV, C Uq and an open set say 
f/jv+ij which is interior to S, i.e. Uj C ^i- Fix a partition of unity {'/^j},-^"'^ subject 

to the covering {UjjfJo^ of n, i.e. ipj G 2?(R"), < ipj < 1, and Ej^ljVi = 1- 
Note that (^q = 1 in a neighborhood of dil. Let (^j denote cut-off functions with 
support in Uj such that ipj = 1 on the support of (pj , and set bj = hipj . 

Let z :— (li, TT, q, h) with g = [tt] be a solution of (13. ip where we assume w.l.o.g. 
uq ^ ho = hi ^ f ^ fd ^ 0, and (ubli^of^) = I(ws|i^s)l = 0. We then set Uj = ipjU, 
TTj = ^jTT, qj = ipjQ, hj — ipjh, as well as u^j = ipjUi,, u^j = (pjU-^, gj = Vjg, and 
ghj = fjgh- Then for j = 1, . . . , N, the quadruples Zj := (uj, ttj, qj, hj) satisfy the 
problems 

pdtUj — iJ^Auj + Vttj = Fj{u,tt) in M"\I]j, 

div Uj ^ (Vfjlu) inM"\i;j, 

|-/i([Vuj] + [Vuj]^) + Qjjiysi - cr{^Sjhj)iyj:. = g^ + Gjiu) on E^-, 

Kl = "Sj, q] = [tTjI on Sj-, 

dthj - {uj\vY.j) + {bj\Vy,hj) = ghj + Ghj{h) on S^-, 

Uj{{)) = in R"\Ej, /ij(0) =0 on S^-. 
Here we used the abbreviations 

Fj{u, tt) = [V<^j]7r - ^[A, ^j>, 

Gj{u) = [-^l{V(pJ ®u + u® Vipj)\v^. - cr[As, V3j]/i!^Sj, 
and 

Ghj{h) = {hj\VT.^j)h. 
For j — Q we have the standard one-phase Stokes problem with parameters p2, fJ-2 
on n with Dirichlet boundary conditions on dVl, i.e. 

P2dtUo ~ ^2Auo -I- VtTj = i^ol"", tt) in fi, 

div Mo = {Vipo\u), in 51, 

Wo = "60 on dfl, 

uo{0) =0 in rj. 

For j = A^ + 1 we obtain the one-phase Stokes problem on R" with parameters 
Pi, Pi, i-c. 

pidtUN+i - HiAuN+i + Vttat+i = FAr+i(u,7r) in R", 

div un^i — (V(y9jv+i|w) in R", 

ujv+i(0) = inR". 

Concentrating on j = 1, . . . , A'^, we first note that [A, ipj] are differential opera- 
tors of order 1, hence if u e qZ^ then 

[A, if,]u G oi?pi/2( J; Lp(R")") n ip(J; i7pi(R"\I],)"). 

Since f = fd = the pressure tt belongs to 

vr e oi?p(J;ip(R")) nip(J;i7pi(R"\S,)), 

by CoroUarv 13.21 hence we have 

F,{u,tt) e oi?p(J;ip(K"))"nLp(j;i?pi(R"\i],)r, 
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for some fixed < a < | — ij-. Similarly we have 

and since [Aj] ■ ,(pj] is of order 1 as well, we obtain 

This shows that we have 

The terms Ghj{h) do not have more regularity, however, the Banach algebra prop- 
erty yields the estimate 

\Ghjih)\YO < C\b\Yo\h\YO < C\b\Yoa^\h\z^, 

with an appropriate exponent 7 > 0. Next we decompose 

Fj {u, n) = Fj (u, n) + Vijjj , 

such that divFj(u,7r) = in W\T,j and (|Fj(u,7r)]|zys,) = on E^. Thus Fj{u,tt) 
is the Helmholtz projection of Fj{u, it) in R". Then 

Also, we decompose Uj = Uj + V^j, where 0^ solves the transmission problem 

A<?!)j = (V(pj» inR"\Ej, 

|p0,]=O onS,, 

ia,,^^,]=0 onS,. 

Note that 

V0, G oi?^(J;i?p'(M"\S,)") nL,(J;if3(K"\s^.)«), (3.8) 

by Theorems 18.11 and 18.61 since Ej is smooth. The jump of its trace on E^ then 
belongs to 

[V</>,1 G oHliJ; W^;-i/^(E,)") n Lp( J; Wl-'/^i-E^T). 

and its normal part vanishes, by construction. Further we have 

[mV^.^,] G oWl-'/'^{J- Lp(E,)"X") n L,{J- t¥p2-i/p(s,)"X"). 

Then we set 

TTj = TTj — Ipj + pdfCpj — AtA(/>j, 

and we observe that on Ej 

since by construction ipj and p0j have no jump across Ej. Now the quadrupel 
Zj := {uj,TXj^qj,hj) satisfies the problem 

pdtiij — p.Auj + VtTj — Fj{u,tt) in M"\Ej, 

div Uj ^ {\7ipj\u) inM"\Ej, 

|-/i([Vu,] + [V7i,]T) + q^p^^ - f7(As,./i,>s, = 5j- + G,(u) on E,, ^ 

[%■] = "sj - lV0j], gj = fc-1 onEj-, 

dthj - {uj\v^j) + (bjlV^hj) = ghj + Ghj{h) on E^-, 

Uj{0) = 0, in R"\Ej, /ij(0) =0 on E^. 



16 MATTHIAS KOHNE, JAN PRUSS, AND MATHIAS WILKE 

Here Gj and Ghj are given by 

G,{u) = G,{u) + 2|AiV20,l^s, - lAi(V(^,»l^.s, 
and 

Ghjih) = Ghj{h) + a^Sj^j. 
For the remaining charts with index j = 0, A'^ + 1, i.e. the one-phase problems, the 
procedure is similar. In case j = we use the regularity 

instead of p.Sp . 

We write p.9p abstractly as 

and by Theorem 13.11 for bent interfaces we obtain an estimate of the form 

I^Ie < Ca{\Hj\r + \Bjz\r), 

with some constant Co independent of j. Here E means the space of solutions and 
F the space of data. Since all components of BjZ (except for Ghj{h)) have some 
extra regularity, there is an exponent 7 > and a constant Ci independent of j 
such that 

\Bjz\r<a^Gi\z\M. 
In addition, by Corollary 13.21 we obtain 

\n,\H^ij,LAn))<C2m\r + \B,z\r)<C2m\r + a^'CiC2\z\E). 

This in turn implies 

\9t(f>j\H^{j;Lp(n)) < C2|-ffj|F + a''C3|z|E, 

and then also 

\zj\e < CilHjlw + o'CqIzIe. 
Summing over all j yields z = '^j Zj, hence 

\z\E<C^\H\E + a^CM^- 
Therefore, choosing the length a of the time interval small enough, we obtain the 
a priori estimate 

\zW<Cs\H\r. (3.10) 

Since the equations under consideration are time invariant, repeating this argument 
finitely many times we may conclude that the operator i : qE ^ oli^ which maps 
solutions to their data is injective and has closed range, i.e. i is a semi-Fredholm 
operator. 

It remains to prove surjectivity of i. For this we employ the continuation method 
for semi-Fredholm operators. The estimates are uniform in the densities pj and the 
viscosities /ij, as long as these parameters are bounded and bounded away from 
zero. Hence L — L{pi, p2, fJ-i, fi2) is surjective, if L(l, 1,1,1) has this property. 
Next we introduce an artificial continuation parameter r e [0, 1] by replacing the 
equation for the free boundary h with 

dth + r(-AE)i/2ft - (1 - T){{u\iy^) - (&|Vs/i)} = gn on E. 

The arguments in [ini HO] show that the corresponding problem is well-posed for 
each T 6 [0, 1] in the case of a flat interface, with bounds independent of t £ [0, 1]. 
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Therefore the same is true for bent interfaces and then by the above estimates also 
for a general geometry. Thus we only need to consider the case pi = P2 = Mi = 
P2=T =1. 

To prove surjectivity in this case, note that the equation for h is decoupled from 
those for u and tt, and it is uniquely solvable in the right regularity class because of 
maximal regularity for the Laplace-Beltrami operator. So we may set now h = 0. 
Next we solve the parabolic transmission problem to remove the jump of u across S 
and the inhomogeneity g in the stress boundary condition. The remaining problem 
is a one-phase Stokes problem on the domain 51, which is well-known to be solvable. 
This shows that we have surjectivity in the case pi = P2 = fJ'i = fJ-2 = t = 1, hence 
also for arbitrary p, fi and t = and the proof of Theorem 13. H is complete. 

We close this section with a remark on the situation, which occurs in the treat- 
ment of two-phase flows with variable surface tension. To be precise, Theorem 13. II 
may be generalized to this situation by means of the following 



Corollary 3.4. The statements of Theorem l 3. 1\ remain valid, if the surface tension 
IS not a constant but a function a € C°^'^/'^{J, BC{J:,m+)) n BC{J,C°^'^{J:,R+)). 

This generalization is possible, since the variability of the surface-tension may be 
handeled by a freezing technique during the localization procedure in the previous 
paragraphs. The only difference compared to the case of a constant surface-tension 
occurs in the jump condition of the normal stress in problem (J3.7I) . which has to be 
modified to 

l-piVuj + [Vuj]^)p^^ + qjv^^ ~ crj(As,/ij>s, = Qj + Gj{u) + G„j{h) on S^ 

with aj = (t(0, Xj) > and 

Now, if the radius r > of the balls Uj, j = 1, ... ,N is sufficiently small, we may 
use the Holder/Lipschitz constants Li, L2 > of cr to estimate 

with exponents 71 , 72 > and a constant Cg > depending only on the spatial 
dimension n and the geometry of E. Since Y^ is a Banach algebra, we obtain 

and the a priori estimate p.lOp remains valid, if both a and r are choosen sufficiently 
small. 

4. Local Well-Posedness 



We now turn to problem (jl.ip , and show its local well-posedness for given initial 

Lta To e T 
conditions 



data To S Wp and uq G Wp (Sl\ro)", which are subject to the compatibility 



div uq = in r2\ro, Uq = on dQ, 

|7'r„M(Vuo + [VuoV)iyroJ = 0, {uoj = on Tq, 

where Vvo — I ^^Tq x ^Tq- According to the considerations in Section[2]we approxi- 
mate To for any prescribed 77 > by a real analytic hypersurface E, in the sense that 
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rfi/(A/'^5],7V^ro) < 1], and Fq is parametrized over S by /lo e Wp^ (S)- Employ- 
ing the transformation from Section [5] to the fixed domain, it is sufficient to prove 
the local well-posedness of the quasi-linear problem (j2.2[) . We keep the notation 
and denote by u and n the transformed velocity field and pressure, respectively. 

We are interested in solutions of (12.21) having maximal regularity, and hence, we 
determine a first approximation of the local solution of (j2.2|) by using Theorem l3.1l 
Since the time traces 

(V|uo) e %i-2/P(f]\S) and V,:MVuo + [Vuo]^)z^e1 G W^-'^p{J:) 

will not be trivial due to the transformation, and since the compatibility conditions 
imposed in Theorem l3 . 1 I have to be satisfied, we must be able to construct extensions 
in the right regularity classes to be used as the right-hand sides in ()3.ip . 



Proposition 4.1. Let p > 3, dfl £ C^ , and set J = [0,a]. Let E C fi 6e a closed 
hypersurface of class C^ . Then /o € H-'^{n) n Wp~^^P{n\^) and go G M^p~^/P(E) 
admit extensions 

f e H'p{J;H^\n\j:))nLp{J:H^{n\^)) and 

with /(O) = /o and g{0) = go. 

Proof. We take (^o e H^{^) n Typ~^/''(r2\I]) to be the unique solution of 

A0O = /o in r2\I], 
|(/.o] =0 on S, 
|5,0o] - on S, 
00 = on 9il, 

which exists due to /o G iJj7i(f7)nT4^p"^/^(f7\E); cf. Theorems |H3] and iH Setting 
vq := V0O G Wp^^'^(rj\E)", the parabolic problem 

dtv- Av ^ in J X f^\E, 

H = on J X E, 

|9,i;l = e^^*|a,i;ol on J x S, 

-y = e'^3"*(-!;olao) on J x 9f7, 

t;(0) — Vq in r2\S], 

where As, respectively on Ago denotes the Laplace-Beltrami operator on S respec- 
tively on, admits a unique solution v € H^{J;Lp{n)") D Lp(J; i/2(f7\E)"). We 
may now define 

/ :=divt- e i/pi(J;iJp-i(r!))nLp(J;Hpi(r!\S)) 

and by construction we have 

/(O) = divi;o = A0o = /o. 
Finally, we define g by means of 5 = e^^^go to the result 

g g |^V2-i/2p(j^ ^^(s)) n Lp(j, %i-i/p(i])), 

by the properties of the analytic semigroup e^^*. D 
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Next we introduce the linear operator L = (Li, . . . , L4) defined by the left-hand 
side of (l2?2]). i.e. 



Li{u, vr) := pdtu ~ fiAu + Vtt, 
L2N :- (V|m), 
L3{u, q, h) := l-n{Vu + [Vuf]]!^^ + qv^ - {Aj:h)iy^, 
Uiu, h) := dth - (m|z/s) + (fo| Vs/i), 
and the nonlinearity N — {Ni, . . . , A'4) defined by the right hand side of 

Ni{u,TT,h) := F{h,u)Vu + Mi{h) : V^u 4- Afi(/i)V7r, 
N2{u,h) := Mi{h) : Vw, 

A^3(", h) := Gr{h)Vu + (G^(/i)Vw + G^{h))v^, 
Ni{u, h) := {[Mo{h) - I]V^h\u) + {u- 6|VsM- 
For J = [0, a] let the solution spaces be defined by 

El (a) :={Meffpi(J;Lp(f])")nLp(J;il2(f7\S)") : u = Q on dVl, H = 0}, 

E2(a):=Lp(J;iJpi(rj\I])), 

E3(a) := l^pi/2-i/2p( J. i^(s)") n L,{J- W^-^'Pi^T), 

E4(a) := W^p2-i/2p(j. ^^(j.)) n ^^i(j. M^^2-i/p(5.)) ^ j^^^j. wl-^/^{Y.)). 

We abbreviate 

E(a) := {(ii,7r,g,/i) G Ei(a) x E2(a) x E3(a) x E4(a) : {tt} = q}, 

and equip Ei(a), £2(0), £3(0) and £4(0) with their natural norms, which turn them 
into Banach spaces; E(a) carries the natural norm of the underlying product space. 
A left subscript always means that the time trace of the function is zero whenever 
it exists. Furthermore, the data spaces are defined by 

Fi(a):=Lp(J;Lp(rj)"), 

F2(a) := Fpi(J;F-i(f7)) nLp(J;i/pi(r!\S)), 

F3(a) := W^^i/2-i/2p(j.i^(s)") nip(J;M^pi-i/P(I])«), 

F4(a) := l¥pi-V2p(j;Lp(S)) nLp(J; W-p2-i/p(5.)) 



and 



F(a) := Fi(a) x ¥2(0) x F3(a) x ¥4(0), 



where we equip these spaces again with their natural norms. The generic elements 
of F(a) are {f,fd,9,9h)- 

To shorten the notation we set z — [u, tt, q, h) G E(a) and reformulate the quasi- 
linear problem (|2.2p as 

Lz = N{z), (7/(0),M0)) = (wo,/^o). (4.2) 

From Section [3] we already know that L : E(a) -^ F(a) is bounded and linear and 
that L : oE(a) — > oIF(a) is an isomorphism, for each a > 0. 

Concerning the nonlinearity iV, the following result has been shown in [201 Propo- 
sition 4.1] for the case where E is a graph over R". 
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Proposition 4.2. Suppose p > n + 2 and b G F4(a)". Then 

N eC'-'iE{a),¥{a)), a > 0. (4.3) 

Let DN{u,T:,q,h) denote the Frechet derivative of N at (u,7r,(j, /i) G E(a). Then 
DN(u,T:,q,h) G C{oE{a),o^{o-)), and for any number oq > there is a positive 
constant AIq = Mo{ao,p) such that 

\DN{u,TT,q,h)\c(„K{a),o¥{a)) 

< Mo[\b - u\BC{J;BC)n¥4{a) + l(",7r,g,/l)|E(a)] 
+ Mq [(I V/l|i3C(J;SCi) + \h\Ei{a) + \u\bC{J;BC)) |w|Ei(a)] 
+ A'fo[P(|V/l|BC(J;BC))|V/l|BC(,/;BC) + Q(I V/i|bc(,/:BC1) , |/l|E4(a)) l^lkiCa)] 

for all (u, TT, g, h) G E(a) and a/Z a G (0, oq] • Here, P and Q are fixed polynomials 
with coefficients equal to one. 



The proof carries over to the general case considered here. The basic ingredients 
are still the polynomial structure of the nonlinearity N w.r.t. u and it, which is the 
same as in [20], and the embeddings 

with a = 1 — (n + 2)/p > 0, which show that ^3(0) and ¥4(0) are Banach algebras. 
The difference lies only in the operators Mj{h) which are more complicated in the 
case of general domains, but analytic in h. 

Now, we are able to establish local well-posedness. 

Theorem 4.3. Fix p > n + 2, let dft G C"^, and suppose 

To G w^-^/p, uo G w^-^/p{n\ror. 

Assume the compatibility conditions 

div Uo = in ^\T(), uq = on dfl, 

iVroK'^uo + [Vuo]^)vro] = 0, N] ^ on Tq, 

where Vpa — I — ^Tq x ^Tq • 

Then there exists a — a(uo,ro) > and a unique classical solution (u,7r,r) of 
il.l\) on (0, a). The set 

T= U Wxr(i) 

tG(0,a) 

is a real analytic manifold, and with 

I3:^{{t,x) G (0,a) x n, x^T{t)}, 

the function {u, tt) : O — ?> R"+"'^ is real analytic. The transformed solution (u, tt, q, h) 
belongs to the space E(a). 

Proof. We consider the transformed problem. 

Step 1. Let /lo G Wp^'^^^{n) and Uq G VFp"^/^(fi\I])" be given, such that the 
compatibility conditions are satisfied, and |/io|Bt/c2(s) ^ V- Let Jo — [0, ao] and 

9* G W;i/2-V2p( j„; L^(S)) n Lp( Jo; Wl-'IPi-il)) 
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be extensions of 

(V|wo) e Hp\^) n W^-^/''{n\j:) and r^MVuo + [Vuo]^)z^e] e w^p~'/^(s), 

which exist due to Proposition 23] Further choose an extension u G Ei(ao) of uq 
and set 5 = M restricted to [0,ao] x S. With these extensions we may solve the 
problem 

Lz* = iOJ*,g*,0), iu*{0),h*iO)) = K,/io), 
since all regularity and compatibility conditions of Theorem 13.11 are satisfied. 

Step 2. Wc rewrite problem (|4.2p as 

Lz^ N{z + z*)-Lz* :^: K{z), z e o^a) 

and observe, that the solution is given as z = L^^K{z), since Theorem 13.11 implies 
that L : oE(a) —5- o]F(a) is an isomorphism with 

l^"\(oE(a).oF(a)) < ^' ae(0,ao], 

where M is independent of a < Oq. Thanks to Proposition 14.21 and due to K{0) = 
N{z*) — Lz*, we may choose a E (0, Qq] and r > sufhciently small such that 

|^(0)|f(„) < ^, l^i^(^)l£(„E(a):oF(a)) < ^> ^ ^ oIE(a), |z|e(,) < r, 

hence 

which ensures, that L^^K{z) : B^ (0) — > B^ (0) is a strict contraction; see also 
[20j . Thus, we may employ the contraction mapping principle to obtain a unique 
solution on the time interval [0, a]. 



Step 3. By Proposition 14.21 the right-hand side N is real analytic, and hence, we 
obtain analyticity of (u, tt, q, h) in space and time by the parameter trick as shown 
in [H Theorem 6.3]; cf. also [11 E]. D 



At the end of this section, we want to mention an extension of Theorem 14.31 to the 
case of time- weighted Lp-spaces. For this purpose, for /i G (1/p, 1] we define Ep(a) 
by means of 

z e E^(a) ^ t^-f'z G E(a), 

and similarly we define E^(a). Thus Ei(a) = E(a), Fi(a) = F(a). Such time-weights 
are useful to relax the regularity of the initial values, but maintaining the regularity 
of the solution for t > (5 > for arbitrary small positive 5. More precisely, we have 
the following result. 

Corollary 4.4. Fix p > n + 2, /i G (^ + 3^, 1), let dfl G C^ , and suppose 

To G W^+f'-^/P, Uo G M^p2M-2/p(^\r^)n 



are subject to the compatibility conditions (|4.ip . 

Then there exists a = a(uo,ro) > and a unique solution {u,Tr,q,h) G E^(a) of 
the transformed problem (j2.2p , which depends continuously on the initial data. 
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This result is proved in the same way as Theorem 14.31 taking into account that the 
restriction /i > ^ + ^^t^ still ensures the embeddings 

which are crucial for Proposition 14.21 to be valid also in the corresponding time- 
weighted spaces. Further, the result about the linear problem remains valid in 
these time-weighted spaces as well. This can be shown in the same way as in the 
case /i = 1, taking into account that the operator d/dt admits an iJ°°-calculus in 
Lp,fj.{J;X)) with angle 7r/2, provided /j, > 1/p and X is UMD-Banach space; cf. 
[T7] . We refrain from giving more details, here; see also [15] . 

5. Semiflow, Energy Functional and Equilibria 

5.1. The Induced Semiflovir. Recall that the closed C^-hyper-surfaces contained 
in ri form a C^-manifold, which we denote by A4H^{n). The metric on AiH^{n) 
is defined by 

d(I]i,S2) ■.= dH{M^^uN'^^2), ^1,^2 eMH^n). 

The charts are the parameterizations over a given real analytic hyper-surface E, 
as described in Section 2, and the tangent space at S consist of the normal vector 
fields on S of class C^. This way AiH^{^) becomes a Banach manifold. 

Let d^{x) denote the signed distance for S as introduced in Section 2. We may 
then define a level function ip^: by means of 

^j:{x) ^ g{ds{x)), a;GM", 

where 

g{s) ^ s(l - x{s/a)) + x(s/a)sgn s, s € M, 

and X denotes the cut-off function defined in Section 2. Then it is easy to see that 
S = (^2^(0), and VipY,{x) = i'i:{x), for each a; g E. Moreover, /x = is an eigenvalue 
of V^(^s(a;) with eigenfunction ly-^ix), the remaining eigenvalues of V'^ipsix) are 
the principal curvatures Hj{x) of E at a; G E. 

Consider the subset MH^{n, r) of MH^ifl) which consists of all F G MH^{n) 
such that F C n satisfies the ball condition with fixed radius r > 0. This implies 
in particular that dist(F, 917) > r and all principal curvatures of F G M.H^{Q.,r) 
are bounded by r. Further, the level functions c^r = 9 o dr are well defined for 
F G MH'^in, r), and form a bounded subset of C^(ii). The map $ : MH'^{fl, r) -^ 
C^{fl) defined by $(F) — ipr is an isomorphism of the metric space A4H'^{fl,r) 
onto <^{MH'^{n,r)) C C^{n). 

Lets-{n-l)/p > 2; for F G MH^in,r), we define F G Wj;{n,r) ii ipr G W^{n). 
In this case the local charts for F can be chosen of class Wp as well. A subset 
A C Wp{fl,r) is said to be (relatively) compact, if $(A) C Wp{fl) is (relatively) 
compact. Finally, we define 

distH/p'(F,E) := |(/2r - <Pi:\w-{n) 

forF,EG A^7?2(n,r). 

As an ambient space for the phase-manifold VM of the two-phase Navier-Stokes 
problem with surface tension we consider the product space C(n)"- x MH^{fl). 
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We define VM as follows. 
VM := {(w, r) e C(n)" X MH^{n) : u e Wp'^'^{G\TY\ T e T^p^^/^, 

divM = in n\T, u = on c»f7, 7'rb(Vu + [\/uY)\ut = on T}. 

The charts for this manifold are obtained by the charts induced by J^H (fJ), fol- 
lowed by a Hanzawa transformation; see Section 2. 
Observe that the compatibility conditions 

divu = in rj\r, u = on 9il, 
Pr[/i(Vu + [Vu]"^)]t'r = 0, M =0 on T, 

as well as regularity are preserved by the solutions. 

Applying Theoreni l4.3l and re-parameterizing repeatedly, we obtain a local semi- 
flow on VM. 

Theorem 5.1. Letp > n-\-2. Then the two-phase Navier-Stokes problem with sur- 
face tension generates a local semiflow on the phase-manifold VM . Each solution 
(u,r) exists on a maximal time interval [0,i*). 

5.2. The Pressure. The pressure does not occur explicitly as a variable in the 
local semiflow, the latter is only formulated in terms of the velocity field u and the 
free boundary F. Actually, at every instant t the pressure tt can be reconstructed 
from the semiflow. In fact, fix any t £ (0,t*) and consider G HL{il). Then we 
have by the divergence theorem 

{u{t)\Wcl)) L,in) = -(divu|(/.)i,(n) - ([(u|i^)l|<^)L,(r) = 0, 
hence also (dtu\V(f)) ~ 0. This implies, multiplying the momentum balance divided 
by p with V(j) in L2(^) 



V^ 



V^ ) = ( -Am -u-Vu 



Vcj) 



L2{Q.) \P 

On the other hand, multiplying the stress boundary condition by vr yields 

on r. Thus TT must satisfy the following problem. 

(V- VM = (^Am-u-Vw V0) , (t)eH^,{n) (5.2 

V P J Loin) \P J L-,(n) 



H = aHr + {{^l{Vu + [VuY)\vr\vr), x&T. 

Theorem 18.51 implies that this problem has a unique solution tt G iJp(51\r). Thus 
the pressure is uniquely defined (up to a constant) by the semifiow and can be 
obtained by solving the transmission problem ([57 



^{u,T):^-\p^'Ml,^.+a\T{t)\. 



5.3. The Energy Functional. Define the energy functional by means of 

2^1^ ■ "IL2(0) 

Then 

5t$(u,r) + 2|/.i/2^|2^(^)=0, 

hence the energy functional is a Ljapunov functional, in fact, even a strict one. We 
have the following result. 
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Proposition 5.2. Let pi,/ii,cr > 6e constants. Then 

(a) The energy equality is valid for smooth solutions. 

(h) The equilibria are zero velocities, constant pressures in the components of the 

phases, the dispersed phase is a union of nonintersecting open balls. 

(c) The energy functional is a strict Ljapunov-functional. 

(d) The critical points of the energy functional for constant phase volumes are 
precisely the equilibria. 

This result is a special case of O Theorem 3.1]. 

Remark 5.3. (i) Let us point out that in equilibrium the dispersed phase consists 
of at most countably many disjoint balls Bf{.{xi). If there are infinitely many of 
them, then i?i ^- as i — > oo, hence the corresponding curvatures Hi = —(n — l)/Ri 
tend to infinity, as well as the pressures inside these balls. This is due to the model 
assumption that there is no phase transition. On the other hand, phase transition 
will occur at very high pressure levels. To avoid this contradiction, in the sequel we 
consider only equilibria in which the dispersed phase consists of only finitely many 
balls. Note also that the free boundary will not be of class C^ if f^i has infinitely 
many components. 

(ii) There is another pathological case which we exclude in the sequel, namely if 
the dispersed phase contains balls touching each other. This can only happen if the 
radii of these balls are equal, otherwise the pressure jump would not be constant on 
r. Physically one would expect such an equilibrium to be unstable, but at present 
we are not able to handle this case. Observe that also in such a situation the free 
boundary F is not a manifold of class C^ . 

(iii) Of course neither (i) or (ii) occurs if we assume that ili is connected, the 
continuous phase enjoys this property anyway. 



6. The Stability Result 
Assuming, for simplicity, that the phases are connected, we denote by 

S := {(0, Sr{xo)) : xoen,R>0, Br{xo) C n} 

the set of equilibria without boundary contact. Note that £ forms a real analytic 
manifold of dimension n + 1. Here n dimensions come from the coordinates of the 
center xq and one from the radius R of the sphere Sr{xo). 

Fix any such equilibrium (0, S) G £. We consider the behaviour of the solutions 
near this steady state. Suppose p > n + 2, let dfl E C^, and consider initial data 
{uo,To)€rM. 

Here we have to use the full linearization of the problem at an equilibrium (0, E) 
i.e. at (m, h) = (0, 0), and for this reason we have to replace As in the linear problem 
dSIDby 



i?2 
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This results in the problem 

(6.1) 



pdtV — /LtAv + \7q 


= pfv 


in f7\I], 


divw 


= fd 


in f7\S, 


- qjiyj: — a{A^h)i'^ 


= 9 


on E, 


M 


= 


on E, 


V 


= 


on dfl, 


dth - (u i^s) 


= 9h 


on E, 


Wo in rj\I], h{0) 


= ho 


on E. 



«(0) 

It is well-known that Aj: is selfadjoint, negative semidefinite on functions with zero 
mean, and has compact resolvent in L2(S). Aq = is an eigenvalue with eigenspace 
of dimension n, spanned by the spherical harmonics of degree one. A_i = {n — l)/R^ 
is also an eigenvalue, its eigenspace is one-dimensional and consists of the constants. 
As a base space for our analysis we use 

Xo = Lp^Anrxw^-'^p{E), 

where the subscript a means solenoidal, and we set 

Xi = i^p(r!\E)" X 1^^3-1/" (E). 
Define a closed linear operator in Xq by means of 

A{v, h) = i-{fi/p)Av + p-'Vq, -{v\iy^)), 
with domain Xi := D{A) C Xi defined by 

D{A)^ {{v, h) eXinXa:v = Oon dQ, {vj = and 
|7'sm(Vw + Vi;T)t.s] = on E}, 

where as before "Ps means the projection onto the tangent space of E. 

Here q € iJp(ri\E) is determined as the solution of the transmission problem 



P Jl2 \P J L2 

Iqj = l^liiVv + [Vv]^)i^^\i^s)l + aA^h on E, 
which is well-defined (up to a constant) by Theorem 18.51 This implies 

-Vq = Ti (-AtA +T2{{MVv + [VvY)\v^\v^) + aA^h). 

Then with z — {v,h) and / = {fv,9h) as well as zq — {vo,ho), system (|6.ip can be 
rewritten as the abstract evolution equation 

i + Az = /, t>0, z(0) = zo, (6.2) 

in Xo, provided fd^g^O; see also [25'. 

Since by Theorem 13.11 problem (j6.ip has maximal Lp-regularity, the abstract 
problem (|6.2|) has maximal Lp-regularity, as well. In particular, —A generates an 
analytic Co-semigroup in Xo; see e.g. [16], Proposition 1.2. In addition we have the 
following result. 
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Proposition 6.1. Let pi,^i > 0,cr > he constants, p G (l,oo), and let Xq, A, 

Xi := D{A) be defined as above. Then the following holds. 

(a) The linear operator —A generates a compact analytic Co-semigroup in Xq which 

has the property of maximal Lp -regularity. 

(h) The spectrum of A consists of countably many eigenvalues with finite algebraic 

multiplicity and is independent of p. 

(c) A has no eigenvalues A with nonnegative real part other than A = 0. 

(d) \ = Q is a semisimple eigenvalue with multiplicity ri + 1 . 

(e) The eigenspace N{A) is isomorphic to the tangent space Tz*£ of £ at the given 
equilibrium z* — (0, S), where T, — Sii{xo). 

(f) The restriction of e^ to R{A) is exponentially stable. 

This result is a special case of [H Theorem 4.1] . It shows that equilibria (0, S) G £ 
are normally stable, hence allows for the use of the generalized principle of linearized 
stability to obtain our main result on stability and convergence. 

The following result concerns the stationary Stokes problem 



(6.3) 



It is needed in the proof of the main result of this section. 

Proposition 6.2. Let p > 3 be fixed, and assume that pi and p,i are positive 
constants for i ~ 1,2, and that cu > is large enough. Then the stationary Stokes 
problem with free boundary (|6.3p admits a unique solution {u,t:, |7r|) with regularity 

u e w^~^^pin\j:y\ n e wl-^/p{n\j:), h e w^^^^/pi^), (6.4) 

if and only if the data {fd,g) satisfy the following regularity conditions: 

(a) fdew^^-^^'\n\j:)nH-\n), 

(b) g^{gr,g.)eW^~'^P{J:r. 

The solution map [{fd^g) ^^ (^7^1 [""l)] ^'S continuous between the corresponding 
spaces. 



puou — pAu + Vtt 


= 


in 


Vl\E 


divu 


= /d 


in 


n\Y. 


-V^MVu + [Vu]'')P^ 


^ 9r 


on 


s, 


ilp{Vu+[Vu]'')P^\iy^) + l7TJ 


= gu 


on 


s, 


M 


= 


on 


s, 


u 


= 


on 


dn. 



The proof of this elliptic problem is similar to that of Theorem 13.11 we refer to 
pi] for the case of a flat interface. 

The main result of this section is the following. 

Theorem 6.3. The equilibrium (0, S) G £ is stable in the sense that for each 
e G (0,60] there exists (5(e) > such that for all initial values {uq,Tq) subject to 

dist^^3-2/p(ro,i;) < (5(e) and \Wa\\.^2-2/p^^^^^^^.^ < (5(e) 
there exists a unique global solution {u{t),T{t)) of the problem, and it satisfies 
dist^^3-2/p(r(t),E) < e and ll"(i)llvi/2-2/P(s^yr(t)) - ^' ^ ^ 0- 



(6.5) 
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Moreover, as t -^ oo the solutions {u{t),T{t)) converges to an equilibrium (0,Soo) 
in the same topology, i.e. 

^lim dist^3-2/p(r(i),I]oo) + ll"Wllvi/^2-2/p(^^j,(^^j = 0. 

The convergence is at exponential rate. 

Proof. 1. To prepare, as in [21] we first parameterize the nonlinear phase manifold 
locally near (0, E) over 

X^ := {{u,h) G [W^-'^^PinX^] X W^-'^^P{E)]nXa : u = on dQ, 
M = 0, T'sKVu + [Vu]T)];/s = on E}. 

In particular, this will show that X^ is isomorphic to the tangent space Tiqj2)'PM . 
For this purpose fix w > and solve for given z := (u, h) G Br '' (0) the problem 

puju — [il\u + Vtt = in r2\E, 

divfZ = Mxih) : Vw in f}\S 

-Psb(Vu+[Vu]"r)]j/s = G,(/i)Vu on S, 

-(|M(Vf2+[V?2]T)lz^s|i^E) + H=G.(/i)Vu + G^(/i) on E, 

\u\ =0 on S, 

u = on 951, 

where u = u + u and h = h, i.e. h = 0. We write this equation in short hand 
notation as Li^u = N{z + z) in X^ = Wp (ri\S) n Hp{fl). It is easily shown 
that N is real analytic and A^'(O) = 0; see Section |4l L^ is invertible by Proposition 
16.21 hence the implicit function theorem yields a unique solution u — 4'{u, h) G X^ 
near 0. (j> is real analytic as well, and satisfies (p'iO) = 0. Then we define 

^iu,h) = {u,h) + {(l3{u,h),0). 

Obviously, $ is real analytic, $'(0) = /, $(B^^(0)) C VM, and $ is injective._ 

Hence it remains to show local surjectivity near 0. So suppose that z := (u, h) G 
VM has sufficiently small norm. Solving the problem 

puju — iiAu + Vtt = in fl\S, 

divu = Mi(h) : Vu in fl\i;, 

-Pe[m(Vu + [Vu]"r)]j/s = G^(/i)Vu on S, 

-(|m(Vu + [Vw]T)]j.s|j^e) + W = G,(/i)VS + G^ih) on E, 

|u] =0 on E, 

u — Q on (?il, 

by means of Proposition 16. 2| (u, /i) := {u — u, h) belongs to X.y and (j>{u, h) — u, 
showing surjectivity of $ near 0. In particular, VM is a real analytic manifold near 
(0,1]). 



(6.6) 
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2. Let {u,Tr,h) be a solution on its maximal time interval [0,t*). In this step we 
decompose u = u + u, tt^tt + tt, h — h + h, where (u, tt , h) solves the problem 

piuu + pdtu — /iAu + Vtt — Fu (u, tt, h) in r2\E, i > 0, 

divu^ Mi{h) -.Vu in n\I], i > 0, 

-PslM-E^(M)]t's = G'r(/i)Vw on E, t > 0, 

+ G^(/i)-G^(/i-/i) on E, i>0, (6.7) 

|i2] = on E, i > 0, 

u = on dn, t > 0, 

uh + dth - (u|i/s) = {Mo{h)VT.h\u) on E, t > 0, 

w(0) = iIq in ri\E, /i(0) = /iQ on E, 

with zq = (uo, ft-o) = (0(uo, /lo)) 0) and E{u) := Vu + [Vu]'''. Writing this problem 
abstractly as L^^z = N(z + z) , by the implicit function theorem we obtain a unique 
solution z = z(z) in the function space E(a) for each a <t^. Then z is determined 
by the problem 

pdtU — iiAu + Vtt = puju in r2\E, f > 0, 

divu = in f7\E, f > 0, 

-V^Myu + [Vu]"^)]j/s = on E, i > 0, 

-{lp{Vu+[Vu]'')Ps\iys) + m-aAsh = G^{h) on E, f > 0, 

|u] = on E, t > 0, ^ ' ' 

u = on af^, t > 0, 

dth — (u|/^s) = itjh on E, i > 0, 

u(0) = uq in fi\E, h{0) = ho on E. 

The last equation can be rewritten abstractly in Xq employing the operator A 
introduced above as 

k + AS = R{S),t>0, S{0) = So, (6.9) 

where 

R{z) = {uj{I - Ti)u{S) ^ T2Gj{z),ujh{zj). 
Note that z is a causal functional of z. 

3. Problem (|6.9p is of the form studied in [32], where the generalized principle of 
linearized stability is proved for abstract parabolic quasilinear problems of the form 
(|6.9p . The only difference is that here a part of R is nonlocal, but causal in time. 
Therefore we only comment on the required modifications in the proof of Theorem 
2.1 in 22J. For this purpose we decompose 

R{S) - RnlociS) + Rioc{z) := {uj{I - Ti)u{z),^h{z)) + (-T2G^(z), 0). 

Observe that by construction, if z is an equilibrium then z = z^ hence z — Q. 
Therefore the equilibria are determined by the equation Az^^ — Riodz*)- Further, 
we have an estimate of the form 

\Rnlociz)\¥(to) < e|M|E(to)' 
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provided zq is small in the norm of X^ . 

Let P'^ denote the projection in Xq onto the kernel N{A) along the range R{A) 
of A and let P^ = I — P'^ the complementary projection onto R{A). As in the 
proof of Theorem 2.1 of ^2l we parameterize the set of equilibria £ near over 
N{A) via a C^-map x i-^- x + *(x) such that *(0) = *'(0) = 0. By x := P^z and 
y :— P'^z — ^(P'^z) we introduce the normal form of the problem which reads as 
follows. 

x = r(x,y), x(0)=xo, (6.10) 

y + ^"y = 5(x,y), y(0) = yo, 

where 

T(x, y) - P-{R{^ + *(x) + y) - i?(x + *(x))], 

5(x, y) = P'\R{x + *(x) + y) - i?(x + *(x))] - *'(x)r(x, y), 

since P'=i?(x + *(x)) = and P''R{x + l'(x)) = A"*(x) = A''[x + *(x)]. The first 
component of P'^z equals zero, since the eigenfunctions of A for eigenvalue have 
vanishing velocity part. This implies |u|Ej^(tg) < |y|E(to); hence the nonlocal part is 
estimated as 

\R-aloa{i)\v(ta) < £|"|Ei(io) < ^|y|E(to)' 

for any to > 0. The local parts of T and R can be estimated by e|y|E(to)' ^^ ^^ 
the proof of Theorem 2.1 in [22]. Taking these observations into account, we may 
proceed as in [35] to prove global existence of z, its stability and convergence to 
another equilibrium, provided zq is small in X^. D 

7. Global Existence and Convergence 

Again we assume for simplicity that the phases are connected. There are basically 
two obstructions against global existence: 

- regularity: the norms of either u{t) or r{t) become unbounded; 

- geometry: the topology of the interface changes, or the interface touches the 
boundary of f2. 

Note that the phase volumes are preserved by the semiflow. 

We say that a solution (u, T) satisfies a uniform ball condition, if there is a 
radius r > such that r([0, t*)) C MH^{^, r). Note that this condition bounds the 
curvature of T{t), and prevents it to touch the outer boundary dG, or to undergo 
topological changes. 

Combining the above results, we obtain the following result on the asymptotic 
behavior of solutions. 

Theorem 7.1. Let p > n + 2. Suppose that (w, F) is a solution of the two-phase 
Navier-Stokes problem with surface tension on the maximal time interval [0,t,). 
Assume the following on [0,t,); 
(i) \u{t)\^2-2/p + \T{t)\^3-2/p < M < cxq; 
(ii) (u, r) satisfies a uniform ball condition. 

Then i, ~ cxd, i.e. the solution exists globally, and it converges in VM to an equilib- 
rium (0, Too) G £ ■ Conversely, if (u, T) is a global solution in VM which converges 
to an equilibrium (woo,roo) (z £ in VM, then (i) and (ii) are valid. 
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Proof. Assume that (i) and (ii) are valid. Then r([0,iO) C Wp~^^^{n,r) is 
bounded, hence relatively compact in Wp ~ '^^^{fl, r). Thus we may cover r([0, t*)) 
by finitely many balls with centers S^ such that dist|y3-2/p-e (r(t), Ej) < 5 for some 
3 = jX*); ^ S [O;**)- Let Jfe = {i £ [0,t*) : j{t) = fc}; using for each k a Hanzawa- 
transformation 0^, we see that the pull backs {u(i, ■) o Q^ : i e J^.} are bounded 
in Wp (J7\Efc), hence relatively compact in Wp ^{^\Y.k). Employing now 

Corollary 14.41 with /i = 1 — e we obtain solutions (u"'^,r^) with initial configurations 
{u{t),T{t)) in the phase manifold on a common time interval say (0,a], and by 
uniqueness we have (u^(a), r^(a)) — {u{t + a),r(i + a)). Continuous dependence 
implies then relative compactness of {(«(•), r(-)) : < i < i*} in VM, in particular 
t^, ~ oo and the orbit (u,r)(]R+) C VM is relatively compact. The energy is a 
strict Ljapunov functional, hence the limit set w(m, F) of a solution is contained in 
the set £ of equilibria. By compactness uj{u,T) C VM is non-empty, hence the 
solution comes close to £. Finally, we apply the convergence result Theorem 16.31 
to complete the sufficiency part of the proof. Necessity follows by a compactness 
argument. D 

8. Appendix: Transmission Problems 

In this section we provide some results, concerning the existence and uniqueness of 
solutions to the transmission problem 

Xq-Aq^f, xe n\r 
Ipq] =9, a; e f, 

{d„^q\ =hi, xeT, 

Sd,,,q + (1 - S)q ^ h2j, xedn, 5e {0, 1}, 

where A > 0, 

pix) --^ PiXnAx)+ P2Xn2(x), x&n\T, 
and pj > 0. To be precise, we will study (j8.ip in different functional analytic 
settings. We begin by stating the result for the 'classical' case, i.e. if the basic 
space is given by Lp(f2). 

Theorem 8.1. Let f7 C M" open, 1 < p < oo, f e Lp{n), g e Wp~^''^{T), 
hi e Wl~^'^{T) and h2,5 € Wp^^^^'^{dn), 5 G {0, 1} he given. Then, for each 
A > 0, there exists a unique solution q G H^{^\r) of (|8.ip and a constant Ci > 
such that 

\q\Hi(n\T) < Ci (^|/|lp(o) + |g|j^2-i/p(p^ + |/Jilvi/pi-i/P(r) + If^^Mw^-'-'^^ian)) ■ 

If in addition J = [0, a], / — f{t,x), f £ H^{J\Lp{Vl)) and g = hi = h2.s = 0, 
then for each A > 0, there exists a unique solution q G Hp{J; Hp{fl\T)), and the 
estimate 

MHl{j:H^^(n\r)) < C2||/||ffi(j;Lp(o)) 
holds with some constant C2 > 0. 

Proof. The first assertion basically follows from |12) . since the Lopatinskii-Shapiro 
condition is satisfied at F and dfl. The second assertion follows from the first one 
by differentiating (|8.ip w.r.t. t and by employing the uniqueness of the solution of 
TI). □ 
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We will also need a result for the case A = 0. To this end, let il C M" be a 
bounded domain, g = hi = h2,5 = and / G Lp{il). Define As by Asq = — Aq, 
with domain 

D{As) - {q G H^in\T) : (pq] = ld,.,qj = on T, 

Sd^,,q + {l- S)q = 0, on dn}, S G {0, f }. 

Since 

DiAs) -> Lp{n), 

the resolvent of As is compact and therefore the spectral set <t{As) consists solely 
of a countably infinite sequence of isolated eigenvalues. In case S — 1 it can be 
readily checked that is a simple eigenvalue of Ai, hence Lp{rt) — N{Ai) (B R{Ai). 
The kernel N{Ai) of Ai is given by N{Ai) = Kip, where 

'^p{x) := xnd^) ^ XOala;), a; G fi\r. 

and R{Ai) = {/ G Lp{fl) : (/|lLp) = 0}. Therefore (|8T|) has a unique solution 
q G Hp{il\r)QKtp, provided (/|lp) = 0. In case of Dirichlet boundary conditions, 
i.e. 5 = 0, it holds that N{Ao) = {0}, hence or each / G Lp{fl), the system (|8.ip 
admits a unique solution q G Hp{fl\T). 

Theorem 8.2. Let il C M" a bounded domain, 1 < p < oo, / G Lp(ri), g = hi = 
h2 — and A = 0. Then the following assertions hold 

^'n\r) of dll 

2/ 



(1) If 5 ~0, then there exists a unique solution q G Hp{fl\T) of ()8.1 



(2) If S = 1 and (/|lp) — 0, then there exists a unique solution q G Hp{Q\T) 



ip. 



//in addition J = [0, a], / = /(i,x) and / G H^{J;Lp{n)) s.t. f{t,-) G i?(A5) for 
a.e. t G J, then q G i/pi(J; H^in\T) Q N{As)). 

8.1. A weak transmission problem. Here we study the (weak) transmission 
problem 

(Vg|V0)i,(o) = (/|V0)i,(o), <^ G i/pMl7), 
IpqI =9, x&T, 

where fi C K" is open and bounded with dfl G C^. We want to show that this 
problem admits a unique solution q G Hp{fl\T), that satisfies the estimate 



|Vg|Lp(o) < C {\f\L^{n-R«) + IslvKpi 



i-i/P(P) 



provided / G Lp{n;W") and g G Wp ^^^(T). We will first treat the case / = 0, 

(8.2) 



g G Wp (r) and consider the problem 



Xiq\cj,) + iVq\V<f>)L,in) =0, </) G Hl,{Q), 



- g, a; G F. 

with A > 0. Theorem 18.11 then yields a strong unique solution q G Hp{Vt\r) of (|8.ip 
with f — hi — h2 = which is also the unique solution of (|8.2p . This follows from 
integration by parts. Our aim is to derive an estimate which is of the form 

\q\H^{n\r) < C\g\^i-i/p^^y 
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which will be done by a localization argument. For this purpose we consider first 
the following auxiliary transmission problem 

Xq-Aq^f, xe M", 

lpqj=g, xgR"-\ (8.3) 

Id.qj ^h, xe M"-\ 

with data / e Lp{W), g G Wp~^'^{W-'^) and h e Wp^^'^{W), which will play 
an important role in the forthcoming localization procedure. Solve the full space 
problem 

Xq~Aq^f, xe M", 

to obtain a unique solution qi = (A — A)^^/ e i/p(R"). provided ReA > 0. In the 
sequel we will always assume that A is real and A > 1. In particular, it follows that 

A^/'|<7i|l,(M") + |VgiU^(K„) < C|/|^-.(K„), (8.4) 

with some constant C > being independent of A > 1, since 

AV2|(A - A)-i/|l,(m.>) < CAV2||(/ _ ^)i/2(^ _ A)-i||g(i^^^^)|/|^-.(^„) 

<C\\{I- A)V2(A - A)-'/'\\^^,^.^,Jf\^-.^^„^ 

and 

|V(A - A)-V|l,(M") < C\\il A)(A - A)-i|b(i^,i^)|/|^-.(«„) 

- '^l/lif-i(B")' 

since the norm 

||(/-A)"(A-A)-"||B(i^^i^), a £{1/2,1}, 

is independent of A > 1, which follows e.g. from functional calculus. The shifted 
function q2 = q — qi should now solve the reduced problem 

Xq2-Aq2^0, xeW\ 

lpq2J^~9, a;eM"-\ (8.5) 

ld,q2J=h, xeR"-\ 

with a modified function g e Wp^^^^{R"~^). Let x = {x',y) e R" x R and define 
L := {X — A„)^/^, where A„ denotes the Laplacian with respect to the first n — 1 
variables x' and with domain D{L) = iJp(M"^^). Let furthermore 

P{x\y) = P2X{y>o}{x\y) + PiX{y<o}{x',y), (a;',y)eM""^ X M. 
We make the following ansatz to find a solution of (|8.5p 



q^iy) ■■={:,,;'' 'r:: (8.6) 



e ^!'a+, y >0, 

where a_, a+ have to be determined. The first transmission condition in (|8.5p yields 
P2a+ — PiO- = ff, whereas the second condition implies —L{aj^ + a_) = /i, hence 



TWO-PHASE NAVIER-STOKES EQUATIONS 33 

a+ + a^ = -L-^h. Observe that g,L-^h £ Wp^'^^''{R'^-^). Therefore we may 
solve this hnear system of equations to the resuh 

a- = ■ {g + p2L-^h) , a+ = ] {g + p2L-^h) - L-^h. (8.7) 

Pi + P2 Pi + P2 

In other words, the solution of (|8.5p (hence of (|8.3p ) is uniquely determined and 

a^.a^^ e Wp (M"^-'^). Since \Le^'"' ■ |lp(r"-1xr^) is an equivalent norm in 

Wp (M"^^) and the corresponding constants are independent of A > 1, we 
obtain first 

^ I'72|lp(R") = ^^''^\L~^Lq2\Lj,{s.'^) < C (^|a+lvVpi-i/''(R"-i) + \°—\wl-^'''{K.^- 

Concerning Vq2 in Lp(R"), we estimate as follows 

|V;E'g2|Lp(R") < C'|Log2|Lp(R") = C\L^)L~ -^92|lp(R") 

< C||Loi"^||j3(LpXp)|i92|Lj,(R") 

with Lq := (J — Aj;/)^/^. Here the norm ||-Zjo-^^^||e(L ,l ) does not depend on 
A > 1, which is a consequence of the fmictional calculus. The estimate for dyq2 in 
Lp(]R") is even simpler, since 

This yields the estimate 

A^/^|g2|Lp(R") + |Vg2Up(R") < C (l5lvy^i-i/P(R.-l) + l-^"^^l Wpi-i/P(R„-l)) ■ 

For each fixed A > 1 the operator L~^ is bounded and linear from Wp (R"~^) to 
Wp"^/P(R"-i), where Wp^'^{W"-^) is the topological dual space of iyp/^(R"^i), 
and 1/p + 1/p' = 1. We want to show that the bound of L~^ is independent of 
A > 1. This can be seen as follows. We have 

\L'~ /i|vy^i(R,.-i) <C\LqL~ /i|ip(R™-i) <C\\LqL~ ||e(Lp,Lp)|^Up(R"-i) 

which holds for all h G Lp(R"^^), since \Lq ■ \[^ (R"-i) is an equivalent norm in 
Wj^(R"^^). On the other hand we have 

\L^ ^Ilp(K"-i) = \LaLQ L^ ft,|ip(R„-i) = \LoL^ L^ ft,|ip(Rr.-i) 

< l|ioi""^lle(Lp,Lp)|icr^^lLp(R"-i) 
<C\\L^L^ lli3(ip,Lp)|/i|H'p-i(R"-i) 

for all h e Wp^{W^~'^), since \Lq^ ■ |lp(R"-i) is an equivalent norm in Wp"^(]R"~^) 
and since L~^ and Lq are commuting operators. Finally we apply the real inter- 
polation method to obtain 

|L" ^lM/pi-i/P(R„-i) < C|/i|^-i/p^jj,.„ij, 

for all h e Wp '^(R"^^), where the constant C > is independent of A > 1. In 
summary we derived the a priori estimate 

A^/^k2|Lp(R") + |V92|lp(R") < C (l5lvVpi-i/P(R„-i) + l^lvVp-^/^CK—i)) ' 
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W, 



-i/p^ 



for the solutfon of (|8.5|) . hence 

< C* (|/|ff-i(R,.) + l5'l^^i-i/p(M„-i) + \h\ 
for the solutfon of (|8.3p . since 

- l5lH'pi-i/P(B—i) +C'I/Ih-1(R")' 

by (|8.4p . Consider now a bounded domain fi C M" with 9ri e C^ and let F C fi 
be a hypersurface such that F G C^ , F n d^ — and such that F divides the set 
Q, into two disjoint regions J7i, i72, where Si^i = F and dVi2 = dil U F. Since J7 is 
compact, we may cover it by a union of finitely many open sets Uk, k — 0, . . . , N 
which are subject to the following conditions 

• OQcUq and C/o n F = 0; 

• UiCfli and C/i n F = 0; 

• Uk nr ^ $, Uk n dn = d) k = 2, ..., N and 

N 

\JUkDT. 

k=2 

For fc > 2, the sets Uk may be balls with a fixed but arbitrarily small radius r > 0. 
Let {</Jfe}feLo t>e a partition of unity, such that supp ifk C Uk and < ipk{x) < 1 for 
all a; G n. Consider the transmission problem 

A? - Ag = 0, a: e f7\F 

Nl=5, :^eF, 

[9,,g] = 0, X G F, 

9^(7 = 0, a: e 9ri, 

where g € Wp (F). Set qk = qfk and gk — g^k- By Theorem 18.11 there exists a 
unique solution q e H'^{Q\V) of ([SJ]) . if e.g. A > 1. Muhiplying (|8Jl) by (/?o yields 

Ago - Ago = -2(Vg|V<y9o) -9A(po, x e O, 

which is an elliptic boundary value problem in Q,. Denote by (_foi Go) the right hand 
side of (|8.10p . By [551 Theorem 3.3.4], there exists a common bounded extension 
operator E from Lp{^) resp. H7^^{^) to Lp(M") resp. iJ_7^(M"). Solve the equation 

Ag^ - Ag^ =£;Fo, x e M". 

The solution is given by gQ = (A — A)^^_EFo and we have the estimate 

A'/'|'Z^|l,(r..) + \^qlK(m < C|i?Fo|^-i(^„) < C|Fo|^-i(^) < C|gU^(o), 

as we have already shown. Note that since Fq e Lp{n), it holds that 

kolff|(K") = l(A - A)-ii?FoU2(K„) < C\Fo\L,in) < C|gUi(o), 

and C > does not depend on A > 1. In particular, the real interpolation method 
yields 

l9olwpi+^(M") ^ Clllw^m, s e [0,1]. 
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The shifted function q^ = qo — q^ solves the problem 

Ml - Ag2 = 0, xe n, 

with some modified function Gg G Wp (dfl). By fP, Theorem 9.2], there exists 
a bounded solution operator Sq : Wp (dfl) -^ Hp{fl) such that q^ = SqG^ and 
there exists a constant C > being independent of A > 1 such that 

This yields 

A^/^kolipCn) + |Vgo|L^(o) < C'(l(V<?|V(y9o)lH-i(o) + ll^'PolH-^in) 

Since ipo is smooth and compactly supported and since v £ C^{dfl), we have 

A'/'ko|L,(^) + |Vgo|L,(o) < C\q\w^^n), (8.12) 

for some s G (1/p, 1), since 

l^lw-^/'iao) - C'l^UpOn) < C\q\w^^n), s e (l/p,l), 
and 

l^i^lolw-'^^idn) - ^\9uql\Lp{dn) < C|golvi/pi+=(n) ^ C'|q'|vK^=(a). 
In a next step we multiply (|8.9p by (/?! to obtain the full space problem 

Xqi- Aqi = -2{\/q\\/ipi)-qAipi, x eW\ (8.13) 

This problem admits a unique solution qi = (A — A)^^Fi, provided A > 1, where 
5'i = (A - A)-i : iJ-i(M") -J> -ffp^(M") is bounded and Fi denotes the right hand 
side of (|8.13p . As before we obtain the estimate 

A'/'ki|L,(R") + |VgiU^(K„) < C|(?|i^(n), (8.14) 

with C > being independent of A > 1. 

We turn now to the charts Uk, k = 2,...,N. Multiplying (|8.9p by ipk, k — 
2, . . . , iV, we obtain the pure transmission problem 

Xqk - Aqk = -2{Vq\Vipk) - qAifik, X G M"\r, 

{pqkj =gk, x G r, (8.15) 

Id^qk] = Iqjdi.ifk, x eT. 

Let xq G r. Then there exists k G {2, . . . , N} such that xq G Uk- After a translation 
and a rotation of coordinates we may assume that a;o = and that the normal v{x()) 
at Xq which points from J7i to ft2 is given by i^ixo) = [0, ... ,0, —1]^. Consider a 
graph 7] G C^(]R"~^) with compact support such that 

{{x', Xn) G Uk C K"-i X M : a;„ ^ ?7(x')} - T n [/fc. 

Note that, since Vx'Vi^) ~ 0: ^'^ ni^'Y choose |Vj:'?7|oo as small as we wish, by 
decreasing the size of the chart Uk- Let q{x',Xn) — v{x' ,Xn—r]{x')), where {x\xn) G 
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C/fc. We define a new coordinate by y = Xn — ri{x'), (x',x„) £ Uk- Then we obtain 

Aq{x', Xn) = Ayv{x', y) - 2dy iV^>v{x', y)\\7^,r]{x')) 

+ d^vix\y)\V^,r^\^-dyv{x',y)A^,r]{x') 
and 



1 



since the normal at a; G C/^ fl F is given by 

Let (Ou)(a;',y) :— q{x' ,y + ri{x')) = v{x' ,y) with inverse (Q^'^v){x' ,Xn+i) = 
v{x',Xn+i — vi^')) = Q{x',Xn+i)- Applying the C^-diffcomorphism 9 to (|8.15p 
and considering the terms on the right hand side of ()8.15p which depend on u as 
given functions {fk,gk,hk) yields the problem 

Xvk -AyVk = F{fk,Vk,Lpk,ri), {x',y) e M", 

[pvk}^G{gk), x'eM"-\ y = 0, (8.16) 

ldyVkl = H{vk,Vk.v), x'eM"-\y = 0. 

which is of the form ()8.3p . Here 

F{fk,Vk,(pk,ri) ■= -^dyiyx'VklVx'V) + 9^Ufc|Vj,'r?p - dyVkA^-rj, 

G{gk) := 6.gfc and 

H{hk,Vk,ipk,v) ■= i_^_\^ , 12 {l^x'Vk]\\'x'V) 
We want to apply (|8.8p to (j8.16p and estimate as follows. 

\dy (V,,ffc|V,,77) |,i,-i(K„) < C|(/ - A,)-i/2a^ (V,,t;fe|V,,77) |i^(K„) 

< C| (Vaj'-yfclVj;'/;) |lp(B") 

< C\Vx''']\Qo\vk\^n(Jgl^„y 

In the same way we obtain 

\dyVk\Vx'V\'^\w-^{R'^) ^ C:\'^x'V\'L\'"k\wi(Rr.y 

whereas 

since rj is smooth. Concerning the terms in the Neumann transmission condition, 
we obtain by trace theory 

< C\Vx'ri\c"(R"-^)\vk\w^R^)i 

where a e (1/p, 1). These estimates show that the right hand side of (|8.16l) may 
be estimated by terms that are either of lower order or of highest order, but the 
higher order terms carry a factor of the form |Va;'77|^, 9 > 0, which becomes small, 
by decreasing the size of the chart Uk ■ Applying perturbation theory it follows that 
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there exists Aq > 1 such that for each chart Uk, k — 2, . . . ,N, the hnear problem 
()8.16|) has a bounded solutfon operator 

provided A > Aq. This in turn yields that Q^^SkQ is the corresponding solution 
operator for problem (|8.15p . i.e. we have 

qk = ie-^Ske){Fk,Gk,Hk), 

for each k = 2,...,N, where {Fk,Gk,Hk) denotes the right hand side of (|8.15p . 
Since is a C^-diffeomorphism, we obtain the estimate 

A^^^kfc|Lp(f2) + |Vgfc|Lp(fi) < C (l.9lvi/pi-i/P(r) + kkp=(f2\r)) , (8.17) 

for some s G (1/p, 1) and for each k — 2, . . . ,N. Here the constant C > does not 
depend on A > Aq, as we have already shown in the investigation of (|8.3p . Let us 
introduce 

\v\>.,w^in) := \M'^^\v\LAn) + \^v\LAn), A > 1, z; G Wl{n\T), 

which is an equivalent norm in Wp{i}\T). This yields 

N 
k=0 

with constants C, M > 0, being independent of A. Since s £ (1/p, 1) we may apply 
interpolation theory to the result 

\q\wf{n) < £\q\w^{n) +C{e)\q\L^(^n) 
<£\q\\.wiin) +C(e)|<7|L^(o) 

since by assumption A > 1. Choosing first e > small enough and then A > 1 
sufficiently large, we finally obtain the estimate 

\q\w}m <C\g\^i-i/p^^^ (8.18) 

for the strong solution q S Wp{^\T) of (|8.9p . Now we want to reduce the regularity 
of g. Fix g 6 Wp~^'''{T). Then there exists a sequence {g„i) C Wp~^'^{T), 
such that grn -> 5 as m ^ oo in Wp~ (T). We denote by g„i G Wp{n\r) 
the corresponding solutions of (18.91) . Then it follows from (|8.18l) that (q™) is a 
Cauchy sequence in Wp{n\T). Therefore the limit limm^oo ^m —'■ <Zoo exists and 
900 G W^p (^\r) is the unique weak solution of (|8.9p for sufficiently large A > 1. 

Lemma 8.3. Let 1 < p < 00, 1/p + l/p' = 1 and let g £ Wp (F) be given. 
Then there exists Aq > 1 such that the problem 

A(<zl0)L.(n) + (Vg|V,^)i,(f,) =0, e ff^M^), 

HI =9, x€T, 

has a unique solution q G H}^{fl\T), provided A > Aq. Moreover, the solution 
q G i7p(ri\F) satisfies the estimate 

\q\Hiin) <C\g\^i-i/p^^y (8.19) 
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In a next step we consider the problem 

A(9|'^)l.(o) + (V9|V(/.)i,(o) = (/|V0)i,(o), cP&H;,,{n), 
Ipqj = 0, a; e r, 

where / € Lp{n;R'^) is given. Observe that the mapping -0/ : HL{il\T) — > M 
defined by 



is hnear and continuous, since 

1^-/(0)1 < l/lLp(O;a")l0lHi,(n), 
hence V'/ G (iJ^, (f7\r))*. With the help of the Dirichlet form 

a : H^{n\r) X H^,{n\r) ^ M, a{q,v) := f \/q- Vvdx, 

Jn 

we define an operator A : Hp{V,\r) -^ (iJi(ri\r)) by means of 

{Aq,v) ;= a{q,v), 

with domain D{A) = {q e H^{fl\r) : {pq] = on T}. Making use of these 
definitions, we may rewrite (|8.20p in the abstract form 

Xq + Aq = ijf, in (i7pi,(r!\r))* . (8.21) 

Since 

H;,{n\r)^{H;,,{n\r))\ 

the resolvent of A is compact and therefore the spectral set (j{A) consists solely of 
a countably infinite sequence of isolated eigenvalues. By a bootstrap argument it is 
easily seen that the corresponding eigenfunctions are smooth. Hence, defining A2 
to be the part of A in i2(^\r) with domain D{A2) = {q e D{A) : Au e i2(^\r)}, 
it follows that <t{A) — <t{A2). Integrating by parts, we obtain 

D{A2) = {qe Hlin\r) : (pq] = 0, (d.^qj = on T, d,q = on OQ} 

and A2q = -Ag in Q,\T. Let A e <j{-A) = a{-A2) and let q G £'(^2) be a 
corresponding eigenfunction. Then q satisfies the problem 

Xq- Aq = o, X e n\r, 

Ipqi =0, a; e r, 

d„q = 0, a; G d^. 
Multiplying (J8.22P -^ by pq and integrating by parts, we obtain by (|8.22P n 3 ^ 

-A / p\qfdx = — / pqAqdx = — pi / qiAqidx — P2 q2Aq2dx 
Jn\r Jn\r Joi J02 

= PilVqila +/52|V(?2|2 + / {dur<l2P2q2- d^-^qipiqi)dT 

= pi\Vqi\l+p2\Vq2\l+ j d,^q2{pq\dT 
-Pl|V<Zi|^+/92|Vg2|2>0, 
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where qj denotes the part of q in Oj. In particular it follows that A is real and 
A < for all A G a {—A) and if A = then qi and (72 arc both equal to a constant 
in ill and Q2, respectively, satisfying the identity piqi — P2Q2- In other words the 
eigenvalue A = is simple and the kernel N{A) = N{A2) is given by 

N{A)=Ktp, lp{x):^xndx) + —XnM, x e n. 

Therefore, spectral theory implies {H^,{n\r))* = N{A) ® R{A) and H^{n\r) = 
N{A) © Y, where y is a closed subspace of _ffp(J7\r). Note that these decomposi- 
tions reduce the linear operator A. It follows that the equation Aq — F has a unique 
solution qeY C H^{n\T) if and only if i^ e R{A), or equivalently (F, Ip) = 0. If 
c S K, then any other solution q £ iJp(51\r) of Aq = F is given hy q = q + cAp and 
we have the estimate 

|V<7|l,(o) <C\F\^^,^^^^^^y . 

Lemma 8.4. Let 1 < p < 00, 1/p+l/p' = 1 and let f E Lp{n-,W") be given. Then 
the problem 

(Vq|V0)i,(^) = (/|V0)i,(n), e ffpM^), 

= 0, X g r. 



ftas a unique solution q £ iJp(57\r), satisfying the estimate 

|Vq|Lp(o) < C'|/|Lp(0;R") ■ 

For the final step, let v £ Hp{il\r) be the unique solution of 

Ao(^^|0)l,(o) + (Vz;lV</>)i,(o) =0, £ i/p^l^l), 

Ipvj = g, a; £ r. 



which is well-defined, thanks to Lemma 18.31 With the help of this solution v, we 
define a functional V'd G (-f^p'(^\r)) by 

■)/;„(0) := / Vu- V0dx. 



By definition it holds that V'ti(lp) = 0. Since also ipf{lp) = for all / £ Lp{Vt\ 
Lemma 1131 yields a unique solution w £ Hp{Vl\r) of 

IpH =0, .T £ r. 

Finally, the sum g := w + w £ iJp(f2\r) is the unique solution of 

{pq\ =9, X er 

and we have the estimate 
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Theorem 8.5. Letl<p<oo, l/p+l/p' = I, f e Lp(f7;M") and g G Wp~^^^{T) 
be given. Then the problem 

(Vg|V(^)i,(o) = (/|V</>)i,(o), e H^^^m, 

IpqI = 9, X er 

has a unique solution u G iJp(f2\r) satisfying the estimate 

\^q\L^(Q.) < Cl (^|/|Lp(n;Rn) + l5lw-pl-l/P(p) j ■ 

// J == [0,a], / = f{t,x), f e H^{J;Lpin;R")), g ^ 0, then q G H^iJ; H^in\T)) 
and 

8.2. Higher regularity in the bulk phases. The next problem we consider, is 
about higher regularity in the bulk phases J7\r. To be precise, we study the elliptic 
transmission problem 

Xq-Aq = f, xe n\T, 

Ipqj^o, xer, 
ld,,qj = 0, x G r, ^''- -"> 

5d^,,q + {1-S)q = 0, xe dn, S G {0, 1}, 

where / G Lp{il) D Wp{n\T), s > 0, is given and A > 1. It is our aim to find a 
unique solution q G Wp+^(ri\r) of (|8.23p . Note that by Theorem l8T] there exists 
a unique solution q G Hp{^\T) of (|8.23p . Moreover, there exists a constant C > 
being independent of A > 1 such that the estimate 

\q\H}in\r) < C|/U^(a) (8.24) 

is valid. Thus, it remains to show that in addition q G W?+*(i7\r), provided 
/ G Lp{il) n Wp{il\r). For this purpose let dil G C^ and cover the compact set 
n by a union of finitely many open sets Uk, k — 0, . . . , N which are subject to the 
following conditions 

• dncUo and [/q n T = 0; 

• C/i C rJi and C/i n r = 0; 

• C/fe n r 7^ 0, C/fc n 9^ = fc = 2, . . . , iV and 

N 



U Uk D r. 



fc=2 

For fc > 2, the sets Uk may be balls with a fixed but arbitrarily small radius 
r > 0. As before, let {'Pk}^=o ^^ ^ partition of unity, such that supp^Sfc C Uk and 
< Lpk{x) < 1 for all X € Cl. Let qk ■= q^k and fk '■— ffk- 
Multiplying (j8.23p by ipo yields the problem 

A(jo - Ago = /o - 2(Vg|V(/3o) - q^Vo, x eVi, 
Sd^,,qo + (1 - 5)qo = Sqd^,,ipo, x e dn, (5 G {0, 1}. 

Since ipo is smooth and q G Hp{il), the right hand side {Fo,Go) in (18.25^ is in 
W^{n) X Wp^"~^^''{dn), at least for s G (0, 1]. It follows from [30l Theorem 5.5.1 
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& Remark 5.5.2/2] that go e W^+^iil), s £ [0, 1] and 

kolw^+^n) ^ C'(l^o|vi/p=(o) + IGol^y^i+^-i/Pjgfj)) < C [\f\w;(n\r) + l/lLp(n)j , 

by (I8.24p . where the constant C > does not depend on A > 1. Muhiplying (j8.23p 
by (^1 we obtain the fuU space problem 

Xqi-Aqi=fi-2{Vq\\7ipi)~qAtpi, a; £ M", (8.26) 

with a right hand side in W^(M"), s G (0,1], which we denote by Fi. Then the 
solution of (|8.26p is given by 

<7i-(A-A)-iFi. 
If a e {0, 1} and Fi e H^{W') then qi e H^+"{W') and 

kilHr°(M") < ^K^ - A)i+"/2giU^(K.) = C\{I A)i+"/2(A - A)-iFiU^(„„) 
= C\{I- A)(A - A)-i(/ - A)"/2FiU^(r„) 

< C||(/ - A)(A - A)-1||b(l„l,)|(/ - A)«/2FiU^(k„) 

< C'Wil - A)(A - A)-i||6(i^,i^)|Fi|^„(R„), 

since |(/ - A)i+"/2 • |lp(r™) is an equivalent norm in ij2+"(M"), a e {0, 1}. Note 
that the term 

||(/-A)(A-A)-i||^(i^^i^) 
is independent of A > 1, which follows e.g. from functional calculus. The real 
interpolation method and (|8.24l) then yield the estimate 

I9ilwp2+=(R") ^ C'|i^i|w-p=(B") < C {\f\w^(^n\r) + I/Ilp(o) j , 

for s G (0, 1], where C > does not depend on A > 1. Next, we multiply (|8.23l) by 
(fik, k G {2, . . . , N}, to obtain the pure transmission problems 

Xqk - Aqk = fk - 2{Vq\Vipk) ~ qA^k, x G M"\r, 

Ipqk} =0, X G r, (8.27) 

{duqkl = {q\d^Vk, a; e r, 

with some function fk G Wp"(]R"\r) n Lp(R"). For each fixed A; G {2, . . . ,iV} we 
may use the transformation described above, to reduce (|8.27p to the problem 

A^A - AV- = F, x' G R""\ J/ G R, 

IpV1=0, x'GR"-\y-0, (8.28) 

ldy^=G, x'GR"-\y-0, 

with given functions F G Wp^(R") and G G wl'^'^^'^iW'-^), s G (0,1]. First we 
remove the inhomogeneity F. To this end we solve the Dirichlet problems 

x' em."--\ y>0, V+(x',0) = 0, 

x' £K"-\ y<0, ^p-{x',0) = 0, 
. Let V=^ e W^+^iR") be defined as 

jip+{x',y), 2/>0, 





AV'~A,,^+-a2^+ = F 


and 






AV^-A.-^--^^^- =F 


wher( 


3 F+ := F\y^o and F" := F j^. 



^^{x',y): 

hA (a;',2/), 2/ < 
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Since ip+{x',0) = ip^ix',0) = and {p^pj = {pjip^ = 0, the shifted function 
ip := ip — ip^ solves the problem 

AV^ - At/^ = 0, x' e W-\ y e R, 

M=0, x'eR"-\y^O, (8.29) 

ldy'ipl=G, x' eW'-\ y = 0, 

where G := G - {dy^Jj^l e VFp+""^''^(K"-i). According to dH]) and dHj]) the 
unique solution of ()8.29p is given by 

^(y)^-^L-^hT~^' ^>"' (8.30) 

^' P1+P2 \p2e^«G, y<0, ' ^ 

where L := {\ — A^'Y^"^ with domain D{L) = i7p(R"^^). Assume for a moment 

that G e Typ~^/^(M"-i). Then it follows from semigroup theory and ([OO)) that 
the solution of (|8.29p satisfies the estimates 

as well as 

IV'Ihkru) < C'|G|^^i-i/P(]g„_i), 

where the constant C > docs not depend on A > 1. This can be seen as in the 
proof of Lemma 18.31 Applying the real interpolation method yields 

for some s £ (0,1] and if G G Wp^'^^^'^iW'-^). We have thus shown that the 
transmission problem (|8.28l) has a unique solution ip £ W?+''(IR") if and only if 
F G W^iK"-) and G G Wp^+'"^/''(R"-i). By perturbation theory, there exists Aq > 
1 such that (jOT)) has a unique solution qk G Wp2+''(IR"\r), s G (0,1], satisfying 
the estimate 

\lk\w^+-(R^\r) < G(^l/felTyp=(R"\r) + l(VglV</jfc)|Wp=(R'.\r) + |gA(^fc|H'»(K"\r) 

provided A > Aq . By the smoothness of (pk and by (|8.24p we obtain the estimate 

lftlwp^+=(R"\r) ^ C" [\f\w^(n\r) + klwi+^n}) ^ ^ (l/lw^Ho\r) + I/Ilp(o) j , 

valid for all /c G {2, . . . , A^} and s G (0, 1]. Since {pk}k=o ^^ ^ partition of unity, we 
obtain 

N 
fe=0 

showing that q G T4^p+''(ri\r), s G (0,1]. It is easy to extend this result to the 
case A G [0, Aq). To this end, let / G Lp{n) D W^{n\T) D R{As), s > 0, where 
As : Hp{n\r) — 7> Lp{n) was defined at the beginning of Section 3. Note that 
R{As) = {/ G Lp{n) : {f\lp)2 = 0} if 5 = 1 and A = and R{As) = Lp{n) if either 
(5 = and A > or (5 = 1 and A > 0. Consider the solution q G H^i^XT) of (|g:^ 
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with A G [0,Ao), which is weh-dGfined thanks to Theorem 18.11 and which satisfies 
the estimate (|8.24|) . Rewriting (|8.23|) ^ as 

Xoq - Aq ^ f + {Xa - X)q, 

we may regard the new right hand side / + (Aq — X)q as a given function, say 
/ e W^{n\T), s e (0, 1]. The above result for ([OSl then yields the estimate 

\l\w^+^n\r) ^ C [\f\w^{n\r) + I/Il^co)) 
< C [jf\w-in\r) + |/|Lp(a)J , 

since 

\Q\w^(n\r) = \q\w-{n) < C\q\H}{n) < C\f\L^{n}, 

by the smoothness of q and by (|8?24l) . If ,s > 1 and / e Lp{n) n Wp'(f7\r), then 
q e H^{n\T), since / e Lp{n)r]H^{fl\T). This additional regularity for q and the 
preceding steps allow us to conclude that q G Wp~^^{Q\T), at least for s e [1, 2]. By 
an obvious argument it follows that q G Wp~^^{il\r) for each fixed s > 0, provided 
/ e Lp{n) n Wp*(f}\r). This yields the following result. 

Theorem 8.6. Let fl C M" be a bounded domain with boundary dil G C^+'*, let 

1 < p < oo, s > and f G Lp(fl) n Wp{Q\T). Then the following assertions hold. 

(1) If S — 1 and A = 0, then there exists a unique solution q G W^"*"" (f2\r) 
Kip of ((05)) . provided that (/|lp) = 0. 

(2) // either 6=1 and A > or S = and A > 0, then there exists a unique 
solution q G W^+^iQXT) of (Ig?^ . 

If in addition J = [0,a], / = f{t,x) and f G H^iJ;Lp{n) nW^+^inXT)) s.t. 
fit, •) G R{As) for a.e.te J, then q G H^iJ; W^+'{n\r) N{As)). 
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